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Dedicated to F. V. Atkinson, one of the pioneers of this subject. 

Abstract. We explicitly determine the high-energy asymptotics for Weyl- 
Titchmarsh matrices associated with general Dirac-type operators on half-lines 
and on R. We also prove new local uniqueness results for Dirac-type operators 
in terms of exponentially small differences of Weyl-Titchmarsh matrices. As 
concrete applications of the asymptotic high-energy expansion we derive a 
trace formula for Dirac operators and use it to prove a Borg-type theorem. 



While the high-energy asymptotics, \z\ oo, of scalar- valued Weyl-Titchmarsh 
functions, m^{z, xq), associated with general half-hne Dirac-type differential ex- 
pressions of the form 



and B a self-adjoint 2x2 matrix with real-valued coefficients, B*^"' e L^{[xo, c])^^^ 
for some n G No(= N U {0}) and all c > xq, received some attention over the past 
two decades as can be inferred, for instance, from [32], [56], [60], [61], [97] (and the 
literature therein) , it may perhaps come as a surprise that the corresponding matrix 
extension of this problem, considering general matrix- valued differential expressions 
of the type 



with /,„ the identity matrix in C™, m e N, and B a self-adjoint 2rn x 2™ matrix 
satisfying B^"^ € L^([a;o, c])^™^^™ for some n € No and all c > xq, apparently, 
received no attention at all. (It should be noted that this observation discounts pa- 
pers in the special scattering theoretic case concerned with short-range coefficients 
^(") g L^([xo, oo); (1 -I- |a;|)da;)^"'^^'", where iterations of Volterra-type integral 
equations yield the asymptotic high-energy expansion of M^{z^xo) as |z| ^ oo to 
any order, cf. Lemma 4.1.) This is not because of a lack of interest in this type 
of problem (we will discuss its relevance below), but simply since it is a nontrivial 
one, which, in many of its aspects, must be regarded as more difficult than the 
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corresponding matrix-valued Schrodinger operator case, which in turn, was only 

very recently settled in [20] . The results proven in this paper show that in leading 
order (and independently of the self-adjoint boundary condition chosen at Xq), 

M+{z,xo) , = ilm + 0{l), (1.3) 
|z|— >oo 

where Ce denotes the open sector in the open upper complex half-plane C+ with 
vertex at zero, symmetry axis along the positive imaginary axis, and opening angle 
£, with < e < tt/2. We are interested in proving the asymptotic expansion 
(1.3) and especially in its higher-order analogs in powers of 1/z, under optimal 
smoothness hypotheses on B. Such results arc then also derived for the 2m x 2m 
analog M{z, x) of M+{z, x) associated with Dirac-type operators on M. 

Our principal motivation in studying this problem stems from our general in- 
terest in operator- valued Herglotz fimctions (cf. [17], [40], [41], [43], [44], [45], [46], 
[51], [113]) and their possible apphcations in the areas of inverse spectral theory 
and completely integrable systems. More precisely, using higher-order asymptotic 
expansions of Af_|-(z,a;), one can prove trace formulas for B{x) and certain higher- 
order differential polynomials in B[x) (similar in spirit to an approach pioneered 
in [48] (see also [37], [39]) in connection with Schrodinger operators). These trace 
formulas, in turn, then can be used to prove various results in inverse spectral the- 
ory for matrix- valued Dirac-type operators D = J-^ — B va. L^(R)^™. For instance, 
using one of the principal results of this paper, Theorem 4.7, and its straightfor- 
ward application to the asymptotic high-energy expansion of the diagonal Green's 
matrix G{z,x,x) = {D — z)~^{x,x) of D, the following matrix- valued analog of a 
classical uniqueness result of Borg [15] for one-dimensional Schrodinger operators 
will be proven in in the context of Dirac-type operators in Section 6. 

Theorem 1.1. Suppose that B is of the normal form B{x) = (f;-^ 

with i?i,i(x) and -61,2(2;) self-adjoint for a.e. a; (E K, and assume that D is re- 
flectionless {e.g., B is periodic and D has uniform spectral multiplicity 2m). In 
addition, suppose that D has spectrum equal to M. Then, 

B{x) = for a.e. xgR. (1.4) 

For related results see, for instance, [1], [2], [23], [36], [47], [53], [55]. Inci- 
dentally, the higher-order differential polynomials in B{x) just alluded to represent 
the Ablowitz-Kaup-Newell-Segur (AKNS) or Zakharov-Shabat (ZS) invariants (i.e., 
densities associated with the AKNS-ZS conservation laws) and hence provide a link 
to infinite-dimensional completely integrable systems (cf., e.g., [7], [18], [27], [28], 
[29], [90], [94], [108], [111], [110], [112], and the references therein), especially, hier- 
archies of matrix- valued (i.e., nonabelian) nonlinear Schrodinger equations. 

Although various aspects of inverse spectral theory for scalar Schrodinger, .Ja- 
cobi, and Dirac-type operators, and more generally, for 2x2 Hamiltonian systems, 
are well- understood by now (cf. the extensive list of references provided in [41]), 
the corresponding theory for such operators and Hamiltonian systems with m x m 
matrix-valued coefficients, m g N, is still in its infancy. A particular inverse spectral 
theory aspect we have in mind is that of determining isospectral sets (manifolds) 
of such systems. It may, perhaps, come as a surprise that determining the isospec- 
tral set of Hamiltonian systems with matrix-valued periodic coefficients is a com- 
pletely open problem. It appears to be no exaggeration to claim that absolutely 
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nothing seems to be known about the corresponding isospectral sets of periodic 
Dirac-type operators in the case m > 2. (More or less the same ignorance applies 
to Schrodinger, Jacobi, and more generally, to periodic 2m x 2m Hamiltonian sys- 
tems with m > 2.) Theorem 1.1 can be viewed as a first (and very modest) step 
toward the construction of isospectral manifolds of certain classes of matrix-valued 
potential coefficients B for Dirac-type operators. 

However, asymptotic high-energy expansions for Weyl-Titchmarsh matrices on 
half- lines and on R, their applications to trace formulas for B{x), and the derivation 
of Borg-type theorems for Dirac operators are not the only topics under consider- 
ation in this paper. We also provide a comprehensive and new treatment of local 
uniqueness theorems for B in terms of exponentially close Weyl-Titchmarsh ma- 
trices. More precisely, in Section 5 we will prove the following result (|| • ||c"ixm 
denotes a matrix norm on C™^"*). 

Theorem 1.2. Fix xq € R and suppose that Bj e L^{[xo,xo + i?])2™x2m 
R > 0, posseses the normal form given in Theorem 1.1 a.e. on (xo-oc), j = 1,2. 
Denote by Mj^^{z,x), x> xq, the unique Weyl-Titchmarsh matrices corresponding 
to the half-line Dirac-type operators in L^([a;o, oo))^™ associated with Bj, j = 1,2 
{fixing some self-adjoint boundary condition at xq). Then, 

if for some a > 0, Bi{x) = B2{x) for a.e. x e {xo,xo -\- a), (1.5) 

one obtains 

||Mi,+ (^,aro)-M2,+ (^,a;o)||c".x,„ = 0(e-2M^)») (1.6) 

\z\ — >oo 

along any ray p+ C C+ with < arg(z) < tt {and for all self-adjoint boundary 
condition at xq). Conversely, if m > 1, assume in addition that Bj e L'^{[xo,xo -\- 
^j^2mx2m^ j _ ^ 2. Moreover, suppose that for all e > 0, 

\\Mi,+{z,xo)-M2,+{^,xo)\\cmxr. ^ = 0(e-2i-W(«-^)), £ = 1,2, (1.7) 

along a ray p+,i C C+ with < a.Tg{z) < n/2 and along a ray C C+ with 
7r/2 < arg(2;) < tt. Then 

Bi{x) = B2{x) for a.e. x G [xa,xa + a]. (1.8) 

We also prove the analog of Theorem 1.2 for the 2m x 2m Weyl-Titchmarsh matrices 
Mj{z,x) associated with Dirac-type operators on M corresponding to Bj, j = 1,2. 

In the context of scalar Schrodinger operators, the analog of Theorem 1.2 was 
first proved by Simon [114]. An alternative proof, applicable to matrix- valued 
Schrodinger operators was presented in [50] (cf. also [41]). More recently, yet an- 
other proof was foimd by Bennewitz [13] (following some ideas in [16]). In fact, 
our proof of Theorem 1.2 is based on that of Bennewitz [13] with additional mod- 
ifications necessary to accomodate Dirac-type operators. These results extend the 
classical (global) uniqueness results due to Borg [16] and Marchenko [91], [92] which 
state that half-line m-functions uniquely determine the corresponding potential co- 
efficient. The Dirac-type results such as Theorem 1.2 appear to be new, even in the 
special case m, = 1. Previous results in the Dirac case focused on global uniqueness 
questions only. We refer to Gasymov and Levitan [34] in the case m = 1 and to 
Lesch and Malamud [81] in the matrix case m € N. 
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Next, we briefly sketch the content of each section. Section 2 provides the nec- 
essary background results on Dirac-typc operators and recalls the basic notions of 
Weyl-Titchmarsh theory for Hamiltonian systems on a half-line as well as on R, as 
developed in detail by Hinton and Shaw in a series of papers [62]-[66] (see also [8], 
[71], [67], [68], [73], [74], [78], [79], [80], [82], [99], [110]). In fact, most of these refer- 
ences deal with more general singular Hamiltonian systems and hence we specialize 
some of this material to the Dirac-type operator case at hand. While our treatment 
of Wcyl-Titchmarsh theory in Section 2 is somewhat detailed, the results presented 
appear to be of vital importance for our asymptotic expansions in Sections 3 and 
4. At any rate, we intended to present this material as concisely as possible. 

Section 3 is devoted to a proof of the leading-order for the asymptotic high-energy 
expansion (1.3) of Mj^{z,x) for the Dirac case. We follow the strategy developed 
in the context of matrix-valued Schrodinger operators in our joint paper [20] by 
appealing to the theory of Riccati equations. By doing so, we follow the lead of 
Atkinson who highlighted the importance of Riccati equations, in this regard, first 
in [9], subsequently in [10], [11] and ultimately in the unpublished manuscript [12] 
in which he obtains the leading order for the asymptotic high-energy expansion of 
Mj^{z,x) for the matrix- valued Schrodinger case. 

Theorems 3.4 and 3.6 contain two characterizations of the Weyl disk (cf. Defi- 
nition 2.7). These characterizations provide an answer in Remark 3.7 to a point 
raised in [20] concerning the nature of the Weyl disk. From these characterizations 
of the Weyl disk, we obtain a realization of M^{z,x) as a differentiable function 
of X which satisfies a certain Riccati equation globally and whose imaginary part 
is strictly positive. We observe, in Remark 3.5, that the totality of Weyl disks, 
Dj^{z,x) (cf. Defintion 2.12), represents the phase space for these solutions. Thus, 
the asymptotic expansion we seek, represents the asymptotic high-energy behavior 
for certain solutions of a given Riccati equation. 

Section 4 develops a systematic higher-order high-energy asymptotic expansion 
of -M_|_(z, x) as |z| oo, combining the leading-order asymptotic result in Section 3 
with matrix-valued extensions of some methods based again on an associated Ric- 
cati equation. More precisely, following a technique in [49] in the scalar Schrodinger 
operator context, we show how to derive the general high-energy asymptotic expan- 
sion of M_|_(z,x) as \z\ ^ oo by combining Atkinson's leading-order term in (1.3) 
and the corresponding asymptotic expansion of M+(^, a;) in the special case where 
B has compact support. Section 5 then contains our new local uniqueness results 
for B(x) in terms of exponentially small differences of Weyl-Titchmarsh matrices 
as indicated in Theorem 1.2. Finally, in Section 6 we derive a new trace formula 
for Dirac-type operators D in L^(M)^'", using appropriate Herglotz representation 
results for the diagonal Green's matrix G{z,x,x) discussed in Section 2. More- 
over, we derive the Borg-type Theorem 1.1 for Dirac operators and close with an 
application to the case of periodic potentials coefficients B. 



2. Weyl-Titchmarsh Matrices for Hamiltonian Systems 

We now turn to the Weyl-Titchmarsh theory for Hamiltonian systems as devel- 
oped by Hinton and Shaw in a series of papers devoted to the spectral theory of 
(singular) Hamiltonian systems [62]-[66] (see also [67], [68], [79], [80], [105], [110], 
[111], [112]). Throughout this paper all matrices will be considered over the field 
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of complex numbers C. The basic assumptions throughout are described in the 
following three hypotheses. 

Hypothesis 2.1. Fix m e N and define the 2m x 2m matrix 




(2.1a) 



where Im denotes the identity matrix in C™^"^. Suppose 

A,-fc,B,-feeLi„,(M)-x™, i,A; = l,2 (2.1b) 

and assume 
for a.e. x G M. 

Li\,^(K) denotes the set of locally integrablc functions on R. With M e C"^™, let 
M* denote the transpose, let M* denote the adjoint or conjugate transpose of the 
matrix M and let M > and M < denote nonnegative and nonpositive matrices 
M (i.e., positive and negative scmidcfinitc matrices). Moreover, let Im(M) = (M — 
M*)/{2i) and Re(M) = (M + M*)/2 denote, respectively, the imaginary and real 
parts of the matrix M. 

Given Hypothesis 2.1, our Hamiltonian system is given by 

J^'{z, x) = {zA{x) + B{x))'P{z, x), zeC (2.2a) 

for a.e. a; e M, where z plays the role of the spectral parameter, and where 

nz,x) = (^^iJ^j) , M^,-) e ^aoc(M)'""^ i = 1,2. (2.2b) 

AC\oc{^) denotes the set of locally absolutely continuous functions on M. The 
parameter r in (2.2b) will be context dependent and range between 1 < r < m. 

For our discussions of the Weyl-Titchmarsh theory for the Hamiltonian system 
(2.2), we introduce the definiteness assumption found in Atkinson [8]. 

Hypothesis 2.2. For all nontrivial solutions ^ of (2.2a) with r = 1 in (2.2b), we 
assume that 

I dx'^{z,x)*A{x)^{z,x)>{) , (2.3) 

J a 

for every interval (a, 6) C M, a < 6. 

A principal example of such a system is the Dirac-type system obtained when 

A{x) = hm, (2.4) 

and the subject of the present paper; another example being the matrix-valued 
Schrodinger system, obtained when 

and the subject of [20]. When (2.5) holds, we note that (2.2a) is equivalent to 

-'^'({z,x) + Q{x)il}i{z,x) = zil)i{z,x), (2.6) 
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Mz,x)=ij[{z,x) (2.7) 

for a.e. a; e K. Hypothesis 2.2 clearly holds in both examples. 

Next, wc introduce a set of matrices that will serve as boundary data for sepa- 
rated boundary conditions. 

Hypothesis 2.3. Let 7 = (71 72) with 7j e C™^'", j = 1,2. We assume that 7 
satisfies the following conditions, 

rank(7) = m, (2.8a) 

and that either 

Ini(727i)<0, or Im(727i) > 0, (2.8b) 

where {2i)~^ jjj* = Im(727i). Given the rank condition in (2.8a), we assume, 
without loss of generality in what follows, the normalization 

77* = Im. (2.8c) 

Remark 2.4. With a G C™^^™, the conditions 

aa* = Im, a J a* = (2.9) 

imply that a satisfies Hypothesis 2.3, and they explicitly read 

aia\ + a2Q.2 = Irm Q.2Q.\ — Q.1Q.2 = 0. (2.10) 

In fact, one also has 

a\ai + 0.2^.2 = Im, a20i\ — a\a2 = Q, (2-11) 

as is clear from 

Q^l I Oi\ —0-2 \ J I Ol\ \ ( Cfl Oi2 



-a2 aij \a2 at J ~ ^™ ~ \a*2 al J {-a2 aj ' ^ ' 

since any left inverse matrix is also a right inverse, and vice versa. Moreover, from 
(2.11) we obtain 

a*aJ+Ja*a = J. (2.13) 

With a e C"*^^"* satisfying (2.9), let ^{z, x, xq, a) be a normalized fundamental 

system of solutions of (2.2) at some xq G M. That is, ^{z,x,xo,a) satisfies (2.2) 
for a.e. a; e M, and 

^{z,xo,xo,a) = {a'' Ja*)=('^l "fV (2.14a) 

\"2 "1 y 

We partition ^'(2;, x, xq, a) as follows, 

^(2;, X, xo, a) = {Q{z, x, xo, a) ^{z, x, xq, a)) (2.14b) 

^ /ei{z,x,xo,a) (j)i{z,x,xo,a)\ (2 14c) 

\92{z,x,xo,a) (f)2{z,x,xo,a)J ' 

where Oj{z, x, xq, a) and (/)j(z, x, xq, a) for j = 1,2 arc m x m, matrices, entire with 
respect to z G C, and normalized according to (2.14a). One can now prove the 
following result. 
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Lemma 2.5. Let @{z, x, xq, a) and ^{z, x, xq, a) be defined in (2.14b) with a and /3 

satisfying Hypothesis 2.3 and, with lm{a2C(l) = 0. Then, for c ^ xq, /3$(z, c, .tq, a) 
is singular if and only if z is an eigenvalue for the regular boundary value problem 
given by (2.2a) on [xq, c] ifo xq and on [c, xq] ifc < xq together with the separated 
boundary conditions 

a^iz, xo) = 0, I3^{z, c) = 0, (2.15) 

where lf'(z,2:) = (^"1(2:, 2^)* ip2{z,xyy with il)j{z,-) G AC{[xo,c]) if c > xq and 
^pj{z,-) G AC{[c,xn]) ifc< xa, j = 1,2. 

Note that the regular boundary value problem described in Lemma 2.5 is self-adjoint 
when Im{P2Pt) = 0. 

In light of Lemma 2.5, it is possible to introduce, under appropriate conditions, 

the m, X m matrix-valued fimction, M{z, c, xq, a, /3), as follows. 

Definition 2.6. Let <d{z,x,xo,a), and ^{z,x,Xo,a) be defined in (2.14b) with 
a and /3 satisfying Hypothesis 2.3 and with Im(a2Q!*) = 0. For c ^ xq, and 
l3^{z, c, Xq, a) nonsingular let 

M{z,c,xo, a, (3) = — [/3$(z, c, xq, a)]""'" [/36(z, c, xq, a)]. (2.16) 

M{z,c,XQ,a, P) is said to be the Weyl-Titchmarsh M -function for the regular 
boundary value problem described in Lemma 2.5. 

The Weyl-Titchmarsh M-function is an m x m matrix-valued function with 
meromorphic entries whose poles correspond to eigenvalues for the regular boundary 
value problem given by (2.2a) and (2.15). Moreover, if M € C™^"', and one defines 

t/(.,.,.o,a)= =*(.,.,.o,a)(y , (2.17) 

with Uj{z,x,XQ,a) G C"*^"*, j = 1,2, then U{z,x,xo,a) will satisfy the boundary 
condition at x = c in (2.15) whenever M = M{z,c,Xo,a,l3). Intimately connected 
with the matrices introduced in Definition 2.6 is the set of m x m complex matrices 
known as the Weyl disk. Several characterization of this set have appeared in 
the literature (sec, e.g., [8], [11], [12], [67], [62], [79], [99]). We now mention two, 
and will introduce two others in Section 3 which we use in the derivation of the 
asymptotic expansions that are the subject of Sections 3 and 4. 

To describe this set, wc first introduce the matrix-valued function Ec{M): With 
c xq, z £ C\IR, and with U{z,c,Xo,a) defined by (2.17) in terms of a matrix 
M e C™^™, let 

Ec{M) = (T{xQ,c,z)U{z,c,xo,a)*{iJ)U{z,c,xo,a), (2.18) 

where 

cris,t,z) = <T{s,t)=a{s,t,i), a(^) = a(l, 0, z), (2.19) 

with s ^ and s, t € R. 

Definition 2.7. Let the following be fixed: Real numbers xq and c ^ xo, an mx2m 
matrix a satisfying (2.9), and z G C\M. ^{z, c, xq, a) will denote the collection of all 

M e C™^™ for which Ec{M) < 0, where Ec{M) is defined in (2.18). V{z,c,xo,a) 
is said to be a Weyl disk. The set of M e C™^"* for which E^iM) = is said to 
be a Weyl circle (even when m > 1). 
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This definition leads to a presentation that is a generalization of the description 
first given by Wcyl [119]; a presentation which is geometric in nature, involves the 
contractive matrices V G c™x™^ g^ch that VV* < Im, and provides the justifica- 
tion for the geometric terms of circle and disk (cf., e.g., [62], [67], [79], [99]). 

The disk has also been characterized in terms of matrices which statisfy Hypoth- 
esis 2.3 and which serve as boundary data for the regular boundary value problem 
described in Lemma 2.5 (cf., e.g., [11], [12]). More precisely, one could have used 
the following alternative definition. 

Definition 2.7A. 'D{z,c,xo,a) denotes the collection of all M G cnx™ obtained 
by the construction given in (2.16) where cj^xq, z £ C/M, where a and /? are the 
mx m matrices defined in Hypothesis 2.3 for which 0-(c, a;o, 2;)Im(/32/3i) > 0, and 
lm{a2al) = 0. 

However, in this paper we take Definition 2.7 as our point of departure. 

We note that the Weyl circle corresponds to the regular boundary value problems 
in Lemma 2.5 with separated, self-adjoint boundary conditions. For convenience 
of the reader, and to achieve a reasonable level of completeness, we reproduce the 
corresponding short proof below. 

Lemma 2.8 ([65], [67], [79]). Let M e C"^™, c ^ Xq, and z G C\IR. Then, 
Ec{M) = if and only if there is a (3 & C™^^™ satisfying (2.9) such that 

() = (3U{z,c,xo,a), (2.20) 

where U{z, c) = U{z, c, xq, a) is defined in (2.17) in terms of M. With (3 so defined, 

M = - [P^z, c, xo, a)]-'[pe{z, c, xo, a)], (2.21) 

that is, M = M{z,c,XQ,a, p). Moreover, (3 may he chosen to satisfy (2.8c), and 
hence Hypothesis 2.3. 

Proof. Let z G C\R, and suppose for a given M G C™^™ that there is a /3 G 
(j-mx2m Satisfies (2.9) and such that (2.20) is satisfied. Given that /3J/3* = 

2iIm(/32/3i) = 0, and given that rank(/3) = rank(J/3*) = m, there is a nonsingular 
C G C"^'" such that U{z, c) = J/3*C. Hence, Ec{M) = ia{c, xq, z)C*(3JI3*C = 0. 

Upon showing that /3$(z,c) = /3$(z, c, xo, a) is nonsingular, (2.21) will then 
follow from (2.20). If f3^{z, c) is singular, then there are nonzero vectors v,w G C™ 
such that (3^{z, c)v = 0, and such that $(2:, c)v = jp*w. Let 'I^j = ^j{z, x), j = 1,2 
denote solutions of (2.2a) with z = zj, j = 1, 2. Then, 

{^*M2y = {z2 - zi)^lA^2. (2.22) 

Using (2.22), and recalling that '^{z,x) is defined in (2.14), we obtain 

2ilm(2;) / dxv*^{z,x)* A{x)^{z, x)v = v*^{z, c)*J^{z, c)v (2.23a) 

= w*pJl3*w = 0. (2.23b) 

Thus, by Hypothesis 2.2, lm(z) = 0. This contradicts the assumption that z G C\M. 

Conversely, if EdM) = for a given M G C™^'", then for z G C\IR let /? = 
{Im M*)'i/{z,c,XQ,a)*J = U{z,c,XQ,a)*J. One observes that (2.20) is satisfied 
and that rank(/3) = m. Moreover, = a{xo,c, z)Ec{M)/2 = Im(/32/3f). If for this 
choice of /?, (2.8c) is not yet satisfied, one introduces A = {(3f3*)~^/'^ (3 and observes 
that = MJ{z,c,XQ,a), Im(A2A|) = (/3/3*)-i/2lm(/?2/3if )(/?/3*)"^/^ and that A 
satisfies all requirements of (2.9). □ 
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Next, we recall a fundamental property associated with matrices in ^{z, c, xq, a). 

Lemma 2.9. If M G 'D{z,c,xo,ct), then 

a{c,Xo,z)Im{M)>0. (2.24) 

Moreover, whenever (3 e C"^^™ satisfies (2.9), 

M{z,c,XQ,a,(3) = M{z,c,xo,a,l3)*. (2.25) 

Proof. Let Sl/j = \I/j{z,x), j = 1,2 denote solutions of (2.2a) with z = zj, j = 1,2. 
Then {^^^2)' = {z2 - ^i)!Z'iM<Z'2 as in (2.22). This implies 

2ilm(^) / dxU{z,x)*A{x)U{z,x) = U{z,x)*JU{z,x)\'' 

= 2ilm(M) + U{z, cyjU{z, c), (2.26) 

with U{z, x) = U{z, X, Xo, a) defined in (2.17). Moreover, by the definition of Ec{M) 
given in (2.18), one obtains 

2a{c,Xo,z)Im{M) (2.27) 
= -E^{M) + 2a{c,xo)\Im{z)\ j ds U {z , s)* A{s)U {z , s) . 

J Xq 

By Hypothesis 2.2 and Definition 2.7, one infers that a{c,Xo, z)lm{M) > 0. To 
prove (2.25), let ^{z,x) = ^{z,x,Xo,a), where 'J' is defined in (2.2). Then, by 
(2.22), 

^{z,x)*J^{z,x) = J, (2.28) 
which impHes J^{z.x){^{z,x)J)* = /2m and hence 

^{z,x)J'^{z,x)* = J. (2.29) 

Thus one concludes 

/3$(z, c)(/3e(z, c))* - pQ{z, c)(/3$(z, c))* = /?J/3* = 0, (2.30) 
from which (2.25) follows immediately by Lemma 2.8. □ 

For c > Xq, the fimction M{z, c, xq, a, /?), defined by (2.16), and satisfying (2.24), 
is said to be a matrix- valued Herglotz function of rank m. Hence, for Im(/32/3i ) = 0, 
poles of M{z,c,xo,a, p), c > xq, are at most of first order, are real, and have 
nonpositive residues. Such functions admit a representation of the form 

M{z,c,XQ,a,P) = Ci{c,xo,a,P) + zC2{c,XQ,a, /3) 

+ J dfl{X,c,xo,a,p) (^j^- , 0x0, (2.31) 

where 6*2(0, xq, a,/?) > and Ci(c, xo,a, /3) are self-adjoint m x m matrices, and 
where fl{X, c, xq, ct, P) is a nondecreasing mx m matrix- valued function such that 



f 

J — c 

0((A, /i],c, a;o,Q!,/?) = limlim— / dva{c,XQ)hn.{M{v + ie,c,xo,a,l3)) . 

(2.32b) 



||dO(A, c, xq, a, /3)||c"'X"' (1 + A^)"^ < 00, (2.32a) 

^l.+s 

A+<5 
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In general, for self-adjoint boundary value problems, 0(A, c, xq, a, /3) is piecewise 
constant with jump discontinuities precisely at the eigenvalues of the boundary 
value problem, and that in the matrix-valued Schrodinger and Dirac-type cases 
C2 = in (2.31) (and later in (2.63) and (2.78)). Analogous statements apply to 
— M(z, c, xo, a, /?) if c < xq. For such problems, we note in the subsequent lemma 
that for fixed /3, varying the boundary data a produces Weyl-Titchmarsh matrices 
M{z, c, Xq, ct, P) related to each other via linear fractional transformations (see also 
[46], [51] for a general approach to such linear fractional transformations). 

Lemma 2.10. Suppose a,/3,7 € C'"^^'" satisfy (2.9). Let = M{z,c,xo,a, P), 
and M-y = M{z,c,xo,j,/3)- Then, 

Ma = [-aJj* + a-f*M^][a-f* + aJj*M^]-\ (2.33) 

Proof. Let Ua{z,x) = U{z,x,Xo,a) and U-y{z,x) = U{z,x,Xo,j) be defined in 
(2.17) with M = Ma and M = M^ respectively. Then, 

^ f3Ua{z,c) =-fU^{z,c). (2.34) 

By the rank condition (2.8a), 

Ua{z,c) = Jp*Ca , U^{z,c) = J0*C^ (2.35) 

for nonsingular Ca, C7 G C"'^™. Thus, by (2.14a), and by the uniqueness of 
solution of (2.2a), there is a nonsingular C G C"*^™ for which 




(2.36) 



By (2.13), 

{a* Ja*)-' = (2-=^^) 
and hence, by (2.36) we see that 

7„ = (a7* + aJ-i*M^)C (2.38a) 

M„ = (-aJ7* + a^*M^)C, (2.38b) 

from which (2.33) immediately follows. □ 
Remark 2.11. From the proof of the previous lemma one infers, in general, that 

U^{z, x) = Ua{z, x){a-i* + aJj*M-y). (2.39) 
Moreover, if ao = {Im 0) and 70 = (0 Im) one observes, in particular, 

M(z,c,a;o,ao,/?) = -M(2,fi,.^;o, 70, /?)"'. (2.40) 

We further note that the sets V{z, c, xq, a) are closed, and convex, (cf., e.g., [65], 
[67], [79], [99]). Moreover, by (2.27) and Hypothesis 2.2, one concludes that E^M) 
is strictly increasing. This fact together with Lemma 2.8 implies that, as a function 

of c, the sets T^iz, c, xq, a) are strictly nesting in the sense that 

I>(z, C2, .To, a) C 15(2;, ci, .To, a) for To < ci < C2 or C2 < ci < xq. (2.41) 

Hence, the intersection of this nested sequence, as c ^ ±00, is nonempty, closed 
and convex. We say that this intersection is a limiting set for the nested sequence. 
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Definition 2.12. Let 'D±{z,xo,a) denote the closed, convex set in the space of 

m X m matrices which is the hmit, as c — > ±00, of the nested collection of sets 
^'(z, c, Xo, a) given in Definition 2.7. 'D±{z,xo,a) is said to be a limiting disk. 
Elements of 'D±{z, xq, a) are denoted by M±{z, xo,a) € C™^™. 



In light of the containment described in (2.41), for xq and z e 

V±{z, xq, a) C V{z, c, xq, a), (2.42) 

with emphasis on strict containment of the disks in (2.42). Moreover, by (2.27), 

M e 'D±{z, Xo, a) precisely when Ec{M) < for all c e (a;o, ±00). (2.43) 

The following Lemma appears to have gone unnoted in the literature. 

Lemma 2.13. Let M e C"^™, c ^ xq, and z e C\R. Then, Ec{M) < if and 
only if there is a P € C™^^™ satisfying the condition 

a{c,xo,z)lm{f52f5l)>0, (2.44) 

and such that (2.20) holds with Uj{z,c) = Uj{z,c,xo,a), j = 1,2, defined in (2.17) 
in terms of M. With l3 so defined, (2.21) holds; that is, M = M{z,c,XQ,a,(3). 
Moreover, (3 maybe chosen to satisfy (2.8c), and hence Hypothesis 2.3. 

Proof Let z € C\M, and for a given M G C™^" suppose that there is a /? € C^^^m 
satisfying (2.44) such that (2.20) holds. The matrices Pj, j = 1, 2, are invertible by 
(2.44), and by (2.20) it follows that 

U{z,c)^(''^l'^Au2{z,c). (2.45) 



By (2.18) and (2.45), one then concludes that 

E,{M) = -2a{c,xo,z)u2{z,cyi3-HYn{P2l3l){l3l)-^U2{z,c), (2.46) 

and hence that Ec{M) < whenever (2.44) holds. 

Upon showing that (3<^{z,c) is nonsingular, (2.21) will follow from (2.20). If 

P^(z, c) is singular, then there is a nonzero vector v G C™ such that P^{z, c)v = 0. 
By the nonsingularity of Pj, j = 1, 2, (l)i{z, c)v = — P^^ P24>2{z , c)v, and as a result, 
(2.23a) yields 



2cr(c, a;o)|Im(2;)| / dxv*^{z,x)* A{x)^{z,x)v 

Jxo 

= -2a{c,Xo,z)v*Mz,cyp^'lm{P2Pl){Pl)-'Mz,c)v, (2.47) 

and hence, a contradiction given (2.44) (cf. (2.3)). 

Conversely, if Ec{M) < for a given M G C™^™, then for z G C\M, Uj{z,c), 
j = 1,2, defined by (2.17), are nonsingular. Indeed, if either ui{z,c) or U2{z,c) are 
singular, then there is a. v G C™, v ^ 0, such that v*Ec{M)v = 0, a contradiction. 
Next, let Pi = I„i and let P2 = —ui{z,c)u2{z,c)~^. Then, for these pj, j = 1,2, 
(2.20) holds. Equation (2.46) now implies that a{c, xq, z)lm{p2pi) > (or c ^ Xq 
and z € C\M. For this choice, P does not satisfy (2.8c). However, one can normalize 
P as described in the proof of Lemma 2.8. □ 

Hence by Lemma 2.13 and (2.43), we see that if M e 'D±{z,xo,ce), then for some 
P e C™^2m satisfying (2.44) 

M±{z,xo,a) = M{z,c,xo,a,p). (2.48) 
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Remark 2.14. To the reader of [20], our study of the high-energy asymptotics of 

the Weyl-Titchmarsh M-function for matrix-Schrodingcr operators, wc offer this 
cautionary note: In [20], D{z,c,Xo,a) represents the set of matrices characterized 
by Lemmas 2.8 and 2.13. However, the homeomorphism that exists between the 
contractive matrices V e C™^™, VV* < Im, and the Wcyl disk, D{z,c,xo,a), 
(cf., [65], [67], [79], [99]) shows that those M e C"^™ characterized in Lemma 2.8 
correspond to the set of unitary matrices while those characterized in Lemma 2.13 
correspond to the contractive matrices for which VV* < Im- Hence, Lemma 2.8 
characterizes part of the boundary while Lemma 2.13 characterizes the interior of 
the Weyl disk as it is defined in Defintion 2.7. As a result, the closure of the set 
consisting of those M <E C™^™ characterized by these two lemmas (i.e., those M 
which correspond to VV* < Im, or to VV* = Im) is the Weyl disk. Thus, for 
deriving high-energy asymptotics for M±{z,Xo,a), it is sufiicient to consider the 
subset of the Weyl disk consisting of those matrices, M E ^mx™^ characterized in 
Lemma 2.8 and Lemma 2.13. This was the approach taken in [20]. 

When T>±{z,Xo,a) is a singleton matrix, the system (2.2a) is said to be in the 
limit point (l.p.) case at ±oo. If 'D±{z, Xq, a) has nonempty interior, then (2.2a) is 
said to be in the limit circle (I.e.) case at ±oo. Indeed, for the case m = 1, the limit 
point case corresponds to a point in C, whereas the limit circle case corresponds to 
'D±{z, Xq, a) being a disk in C. 

These apparent geometric properties for the disk correspond to analytic prop- 
erties for the solutions of the Hamiltonian system (2.2a). To recall this correspon- 
dence, we introduce the following spaces in which we assume that — oo<a<6<oo. 



L\{{a,b)) = S^(f) : {a,b) C^™ ^ dx {^{x) , A(f){x))c2m < ooj, (2.49a) 

N{z, oo) = {(j)e L\{{c, oo)) I J(j)' = {zA + B)(t> a.e. on (c, oo)}, (2.49b) 
N{z, -oo) = {(!)& L\{{-oo, c)) I Jcp' = {zA + B)(j) a.e. on (-oo, c)}, (2.49c) 

for some c e M and 2; e C. (Here {cf), V')c" = 't'j^j denotes the standard scalar 
product in C", abbreviating x G C" by x = (xi, • • • , XnY ■) Both dimensions of the 
spaces in (2.49b) and (2.49e), dimc(iV(z, 00)) and Am\c{N{z,—Qo)), are constant 
for 2! e C± = {C e C I ±Im(C) > 0} (see, e.g., [8], [74]). One then observes that 
the Hamiltonian system (2.2a) is in the limit point case at ±00 whenever 

dimc(iV(2;, ±00)) = m for all z G C\M (2.50) 

and in the limit circle case at ±00 whenever 

dimc(A''(.2, ±00)) = 2m for all z G C. (2.51) 

Next we show that the Dirac-type systems considered in this paper are always in 
the limit point case at ±00. Results of this type, under varying sets of assumptions 
on B{x), are well-known to experts in the field. For instance, in the case m = 1 and 
with Bi,2(x) = B2,i{x) this fact can be found in [118]. For B G C(R)2'"x2m ^ 
more general constant matrix A, this result is proven in [81] (their proof, however, 
extends to the current B e Ll^^{M.) case). More generally, multi-dimensional Dirac 
operators with L^Q^(R")-type coefficients (and additional conditions) can be found 
in [83] . A short proof in the case to = 1 has recently been sent to us by Don Hinton 
[58] . For convenience of the reader we present its elementary generalization to to e 
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N below (see also [19] for a sketch of such a proof). After completion of this paper 

wc became aware of a recent preprint by Lesch and Malamud [82] which provides a 
thorough study of self-adjointness questions for more general Hamiltonian systems 
than those studied in this paper. 

Lemma 2.15. The limit point case holds for Dirac-type systems {i.e., for A = 
in (2.2a)) at ±oo. 

Proof. Let {yi{z, x)}i=i,,,,^k and {wj{z, x)}j=i k' denote bases for N{z, ±oo) and 

Af(z, ±oo), respectively. By Theorem 9.11.1 of Atkinson [8], one has k,k' >m for 
z e C\]R. We now assume that k > m. 

One observes that {yi{z, c), . . . , yk{z, c)} and wi(z, c), . . . , Wk' ("z, c)} are linearly 
independent in C^™+^, where k + k' > 2m + 1. Consequently, there is some s G 
{1, . . . , fc} and some re {1, . . . , fc'} such that 

Wr{z,cyjys{z,c)^0. (2.52) 

By Lagrange's identity, 

Wr{z,xyjys{z,x) = Wr{z,c)*Jys{z,c) (2.53) 

is constant with respect to x. On the other hand, an application of Cauchy's 
inequality shows that the left-hand side of (2.53) is in L^((c, ±oo)). By (2.52) one 
obtains a contradiction and hence concludes that 

dimc(iV(^;, ±oo)) = m. (2.54) 

The analogous argument then also yields 

dimc(iV(^, ±oo)) = m (2.55) 

and hence the limit point property of Dirac-type systems with A{x) = hm in 
(2.2a). □ 

Returning to the general case (2.2a), in either the limit point or limit circle 
cases, M±(z, xo, a) G d'D±{z, xo, a) is said to be a half-line Weyl-Titchmarsh ma- 
trix. Each such matrix is associated with the construction of a self-adjoint operator 
acting on L\{[xo, ±oo))nAC([a;o, ioo))^"* for the Hamiltonian system (2.2a). How- 
ever, for those intermediate cases where m < dimc{N{z,±oo)) < 2m, Hinton and 
Schneider have noted that not every element of d'D±{z,Xo,a) is a half-line Wcyl- 
Titchmarsh matrix, and have characterized those elements of the boundary that 
are (cf. [67], [68]). 

For convenience of the reader we summarize some of the principal results on 
half-line Weyl-Titchmarsh matrices next. 

Theorem 2.16 ([3], [17], [51], [62], [63], [66], [78]). Suppose Hypotheses 2.1 and 
2.2. Let z € C\IR, xq € R, and denote by a,-/ € C'"^^™ matrices satisfying (2.9). 
Then, 

(i) ±M±{z,xo,a) is an m X m matrix-valued Herglotz function of maximal rank. 
In particular, 

Im(±M±(^, XQ, a)) > 0, z € C+, (2.56) 
M±{z,xo,a) = M±{z,XQ,ay, (2.57) 
rank(M±(2;,a;o,Q!)) = m, (2.58) 
limM±(A -|- ie, xq, a) exists for a.e. A e M, (2.59) 
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M±{z,xo,a) = [-aJ^* +a'y*M±{z,xo,7)] x 

x[aj* +aJ^*M±{z,xo,^)]~'^. ''^'^^^ 

Local singularities of ±M±{z,xo,a) and ^M±{z,Xo,a)~^ are necessarily real and 
at most of first order in the sense that 

=Flim(ieM±(A + ie,a;o,a)) > 0, A € M, (2.61) 



± lim — — : ^ > 0, A e R. (2.62) 

(ii) ±M±{z, Xq, a) admit the representations 
±M±{z,xo,a)=F±{xo,a)+ [ dfl±{X,xo,a) {{X- z)-' - X{1 + X"")-') (2.63) 

= exp^C±(a;o,a) + J dXE±{X,xo,a){{X - z)-^ - X{1 + X^^)^, (2-64) 

where 

F±(xo,a) =F±(xo,a)*, ^ ||dr!±(A, xq, a)|ic-x™ (1 + A^)"! < oo, (2.65) 

C±{xo,a) = C±{xo,a)*, < S±( • , xq, a) < a.e. (2.66) 
Moreover, 

^±((A, /u],a;o,a:) = limlim— / di'Im{±M±{v' + ie,xo,a)), (2.67) 

S± (A, xo, a) = lim 7r"^Im(ln(±M±(A + ie,a;o, a))) /or a.e. A e R. (2.68) 

eiO 

{Hi) Define the 2m x m matrices 

U±(z,x,xo,a)= (^±A^^^^^o,ay\ = ^(^ a;,a;o,a) f ,^ 
^ ^ \u±,2{z,x,xo,a)J V > > u, J yM±{z,xo,a)J 



di{z,x,XQ,a) <j)i{z,x,Xo,a) 
92{z,x,xo,a) (j)2{z,x,xo,a) J \M±{z,xo,a.) 



(2.69) 



with dj{z,x,xo,a), and (j)j{z,x,XQ,a), j = 1,2, defined by (2.14c). Then, 

/■iboo 

Im(M±(2;,a;o,Q:)) = Ini(2) / dsU±{z,s,xo,a)*A{s)U±{z,s,xo,a). (2.70) 
In the Dirac-type context, where A = /2m, the m columns of U±{z, •, Xq, a) span 

Ar(^,±oo). 

Up to this point, we focused exclusively on Hamiltonian systems and neglected 
the notion of a linear operator associated with (2.2). We did this on purpose as 
the formalism presented thus far is widely applicable and goes beyond the prime 
candidates such as Schrodinger and Dirac-type systems. However, in the remainder 
of this section and for the bulk of the material from Section 3 on, we will focus on 
the Dirac-type case. Thus, in addition to Hypotheses 2.1-2.3, which are assumed 
throughout this paper, we introduce the following hypothesis taylored to these 
occasions. 

Hypothesis 2.17. Assume Hypotheses 2.1 and 2.3 as well as the Dirac-type as- 
sumption (2.4). 



M-FUNCTION ASYMPTOTICS AND BORG-TYPE THEOREMS 



15 



Assuming the Dirac-type Hypothesis 2.17, we now describe the associated Dirac- 
type operator D on R by first introducing the Green's matrix associated with (2.2) 
and (2.4). Define the 2m x 2m matrix G by 

G{z,x,x') = Uzf{z,x,xo,ao)[M_{z,xo,ao) - M+{z,xo,ao)]^^U±(z, a;', xq, ao)*, 

ao = {Im 0), x^x',ze C\R (2.71) 

Next, let (j) e L2(K)2m a^d consider 

Ji}'{z, x) = {zhm + B{x))i){z, x) + (t>{x), z e C\M (2.72) 

for a.e. x G M.. Then, as inferred from [62], [64], (2.72) has a unique solution 
Tpiz, ■ ) e L2(]R)2m PI ACioc(IR)^™ given by 

il>{z,x)= [ dx' G{z,x,x')^{x'). (2.73) 

The Dirac-type operator D in L^(IR)^'" associated with the Hamiltonian system 
(2.2) and (2.4) is then defined by 

{{D-zyi)){x)= f dx'G{z,x,x')ip{x'), ^ e L2(M)2m^ ^ g c\M. (2.74) 

Explicitly, one obtains 

D = J-^-B, (2.75) 

ax 

dom(D) = {cf>e L2(R)2™ I ^ e ACioc(K)'"; (Jcl)' - B^) e L^Rf""}, 

taking into account the limit point property of Dirac-type systems as described in 
Lemma 2.15. Thus, Z) is a self-adjoint operator in L^(]R)^'". 

As described in [62]~[66], the 2to x 2m Weyl-Titchmarsh matrix M{z,Xo,ao) 
associated with D is then defined by 

M{z, Xo,ao) = (Af,jv(z,a;o,ao))^._^,^^ 2 

= [G{z,xo,xo + 0) + Giz,xo,xo-0)]/2, ^ e C\M. (2.76) 

Actually, one can replace ao = {Im 0) by an arbitrary matrix a = [ai a2] € C™^^™ 
satisfying (2.9) and hence introduces 

M{z,xo,a) = {Mjj,{z,xo,a)). .,^^ .^, ^; e C\M, (2.77a) 

Mi^i{z,xo,a) = [M_{z,xo,a) - M+{z,xo,a)]-^ , (2.77b) 

Mi,2(2, xo, a) = 2-i[M_(z, xq, a) - M+{z, xq, a)]-^[M_{z, xq, a) + M+{z, xq, a)], 

M2,iiz, Xo, a) = 2-i[M_(2, Xq, a) + M+{z, xq, a)] [M_(z, xq, a) - M+{z, xq, a)]-\ 

M2^2{z,xo,a) = M±{z,Xo,a)[M_{z,Xo,a) - M+{z,xo,a)]~^Mzf:{z,Xo,a). 

The basic results on M(z, xq, a) then read as follows. 

Theorem 2.18 ([51], [62], [63], [66], [78]). Assume Hypothesis 2.17 and suppose 

that z e C\K, Xq e R, and that a e C™^^™ satisfies (2.9). Then, 

(i) M{z,xo,ct) is a matrix-valued Herglotz function of rank 2m with representation 

M{z,xo,a) = F{xo,a)+ [ dn{X,xo,a) {{X- z)'^ - X{1 + X"^)-^), (2.78) 
= exp (c{xo, a) + J dX T(A, xo,a){{X- z)"^ - A(l + XY^)^, (2-79) 
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where 

F{xo,a) ^ F{xo,a)*, / ||dl7(A, Xq, a)||c2™x2™ (1 + A^)"^ < oo, (2.80) 
Jr 

C{xo,a) = C{xo,a)*, < T{ - ,xo,a) < hm a.e. (2.81) 
Moreover, 

f2((A, /x],a:o,Q;) = limlim— / dvlm{M{u + ie,Xo,a)), (2.82) 
T(A, xo, a) = lim7r"^Im(ln(M(A + is, xq, a))) for a.e. XgR. (2.83) 

{ii) z e C\spec(Z)) if and only if M{z,XQ,a) is holomorphic near z. 

Here spec(T) abbreviates the spectrum of a linear operator T. 

Next, we explicitly discuss the elementary Dirac-type example where A = 72m 

and B = 0. 

Example 2.19. Suppose A = hm, S and let a e C^^^m satisfy (2.9). Then, 

Q(z X X a)= (^ii^'^^^o,a)\ ^ ^alcos{z{x - Xq)) + a^sm{z{x - xo))\ 
y ' ' 0, ) \$2{z,x,xo,a) J cos{z{x — Xo)) — al sm{z{x — Xo)) J ' 

zeC, (2.84) 

^{z X Xo a) = I _ f-a^ cos{z{x - xq)) + a* sin(z(a; - xq))'^ 



4>2{z,x,Xo, a) J \ al cos{z{x — Xo)) + a2 sin(z(.T - xq)) 

zeC, (2.85) 

U,iz,x,xo,a) = fe^iS!'!'!-"!) = f exp(±..(.-xo)), 



u±^2{z, X, xo,a)J \±i{al T ia2)^ 

z e C+, (2.86) 

M±{z, x, a) = ±ilrn, zeC+. (2.87) 

Compared to the case of Schrodinger operators, it is remarkable that M± {z, x, a) 
in (2.87) is, in fact, independent of a. Put differently, in Dirac-type situations, 
M±{z,x,a) may contain no information on the boundary condition indexed by 

In Sections 4 and 5 we will also refer to half- line Dirac operators D^{a) in 
i^([a:o, oo))^™ associated with a self-adjoint boundary condition at xq indexed by 
a € C™^^"* satisfying (2.9), and hence briefly introduce 

D+{a) = J^-B, (2.88) 

dom(£)+(a)) = {(!)& L^{[xo, oo))^™ | (j) € AC{[xo, i?])^" for all R>0; 

a(l){xo) = 0; {Jcj)' - Bcfy e L''{[xo,oo))'^"'}, 

taking into account the limit point property of Dirac-type systems at -l-oo as de- 
scribed in Lemma 2.15. Thus, Dj^(a) is a self-adjoint operator in L^([a;o, oo))^"*. 
In complete analogy one introduces D-{a) in L^((— oo,a;o])^™- 

Next, wc recall a few formulas in connection with Lagrange's identity needed in 
the proof of Theorem 5.3 assuming a G C"x2m satisfies (2.9). Then, explicitly, 
(2.28) and (2.29) read 

6'2{z,x,xo,a)*6i{z,x,xo,a) — 6i{z,x,xo,Oi)*6'2{z,x,xo,a) = 0, (2.89) 
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(f)2{z,x,XQ,a)*(l)i{z,x,xo,a) - (j)i{z,x,xo,a)*<j)2{z,x,XQ,a) = 0, (2.90) 
(j)2{z,x,xo,a)*6i{z,x,xo,a) - (f)i{z,x,XQ,ay62{z,x,xo,a) = Im, (2-91) 
6i{z,x,xo,a)*(f)2{z,x,xo,a) - 92{z,x,xo,a)*(j)i{z,x,XQ,a) = Im, (2.92) 



and 



(j)i{z, X, xq, a)0i{z, x,xo,ce)* — 9i{z,x,xo, aQ)(j)i{z, x,xo,a)* = 0, (2.93) 
4>2{z, X, Xo, a)62{z, X, xo,a)* - 02{z, x, xq, a)(f)2{z, x, xq, a)* = 0, (2.94) 

(f>2{z, X, Xq, a)6i{z, X,XQ,a)* — 02{z, X, Xq, Q;)</)i (z, X, Xq, a)* = Irrn (2.95) 

9i{z,x,xo,a)(l)2(z,x,xo,a)* - (jjiiz, x, xq, a)92{z, x, xq, a)* = Im- (2.96) 

Finally, we note the connection between $ defined in (2.14b), for different boundary 
value data a, 7 e C"^^™ satisfying (2.9), namely 

^{z, X, Xq, 7) = ^{z, X, Xq, a)a'y* + 0{z, X, Xq, a)aJ^* . (2.97) 

This connection formula follows by the uniqueness of solutions of (2.2) and by the 
identity given in (2.13). It is needed in the proof of Theorem 5.3. 

3. The Leading Order Term in the Asymptotic 

Expansion of M±{z,x,a) 

Assuming Hypothesis 2.17, the principal result proven in this section will be the 
following leading-order asymptotic result for half-line Weyl-Titchmarsh matrices 
M±{z,xo,ao) associated with the Dirac-type operator (2.75), 

M±{z,xo,ao) = ±iIm + o{l). (3.1) 

\z\ — ^00 

Here ao = {Im 0) G £mx2m^ (j^ ^ denotes the open sector with vertex at 
zero, symmetry axis along the positive imaginary axis, and opening angle e, with 

< e < it/2. 

This particular topic originates with the order result of Hille [57] and the as- 
ymptotic formulas of Everitt [30] and of Everitt and Halvorsen [31]. By appealing 
to the theory of Riccati equations, Atkinson in [9], [10], and [11] obtains results 
like those of Hille, Everitt, and Halvorsen, both for the Schrodinger case as well 
as for the scalar-Dirac (m = 1) case. Through a deeper understanding of the role 
played by Riccati theory, Atkinson obtains the first order asymptotic expansion of 
M+(z, X, ao) for the matrix-valued Schrodinger case in an unpublished manuscript 
[12]. Our strategy of proof for (3.1) is patterned after Atkinson's approach which 
also appears in our recent work on the full asymptotic expansion for M_|_(2, a;, ao) 
in the matrix- valued Schrodinger case [20]. 

We begin our discussion by noting two additional characterizations for the Weyl 
disk, T>{z,c,XQ,a), for the general Hamiltonian system (2.2a). 

Lemma 3.1. Assume Hypotheses 2.1 and 2.2. Let z G C\R, c ^ xq, and define 
U{z,x,XQ,a), m terms of M € £jnxm (-2.17). Then M e V{z,c,xo,a) if and 
only if 

a{c, Xo, z)lm{ui{z, X, Xo, a)*U2{z, X, Xo, a)) > 0, x G [xo,c), (3.2) 
or equivalently, if and only if 

a{c,xo, z)lm{u2{z,x,xo,a)ui{z,x,xo,a)~^) > 0, a;e[a;o,c). (3.3) 
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Moreover, M G T>±{z,xo,a) if and only if (3.2) and (3.3) hold for c = ±oo. 

Proof. Let U{z,x) = U{z,x,xo,a), and let Uj{z,x) — Uj{z,x,Xo,a), j = 1,2 with 
xe [xo,c). By (2.26), 

2o-(c,a;o)|Im(0)| / dsU{z,s)*A{s)U{z,s) 

J X 

= (T{xo,c,z)U{z,sy{iJ)U{z,s) ^ (3.4) 

X 

By (2.18), this yields 

2a{c, xo, z)Im{ui{z, x)*U2{z, x)) 

= -Ec{M) + 2a{c,xo)\Im{z)\ f dsU{z, s)*A{s)U{z, s). (3.5) 

J X 

The integral expression in (3.5) is strictly positive by Hypothesis 2.2. This yields 

the equivalence of —Ec{M) > 0, and hence of M G c, xq, a), with the condition 
given in (3.2). The equivalence of (3.2) and (3.3) follows from the observation that 

lm{u2{z, x)ui{z, x)^^) = {ui{z, x)^^)*lm{ui{z, x)*U2{z, x))ui{z, x)^^ . (3.6) 

The analogous characterization of T>± (z, xq, a) now follows from Definition 2.12. □ 

In Lemma 3.1, Mj(z, c), j = 1, 2, are well-defined and Ec{M) = precisely when 
a{c,xo,z)lm{ui{z,c)*U2{z,c)) = 0. A similar statement might not hold for (3.3) 
since ui{z, c, Xq, a) might be singular. In part, the latter point motivates the next 

characterization of the disk. 

Lemma 3.2. Assume Hypotheses 2.1 and 2.2. Let z G C\M, c ^ xq, and define 
Uj{z,x) — Uj{z,x,xo,a), j = 1,2, by (2.17). Then M e 'D{z,c,xo,a) if and only if 

ui{z,x) — ia{c,xo, z)u2{z,x) (3.7) 

is nonsingular for x £ [xq, c] and 

'&{z, x) = '&{z, X, Xq, a) = [ui{z, x) + ia{c, xq, z)u2{z, x)] X 

X [ui{z,x) — ia{c,xo, z)u2{z,x)]~^ 

satisfies 

lm-'&{z,x)*'&{z,x)>0, x€[xo,c), (3.9) 

with nonnegativity holding at x = c. Moreover, M G 'D±{z,Xo,ct) if and only if 
(3.8) is well-defined on [xo,±oo) and (3.9) holds for c = ±oo. 

Proof. Let M e 'D{z, c, xq, a) and suppose that Ui{z, ^)v = ia{c, xo, z)u2{z, for 
^ e [xo,c] and vgC'^,v ^0. Then, 

v*a{c,xo,z)lm{ui{z,£,)*U2{z,0)v = -v*ui{z,^)*ui{z,^)v. (3.10) 

By (3.2), an immediate contradiction results if ^ ^ c. However, if ^ = c, then 
either v*E^{M)v > or Uj{z,c)v = 0, j = 1,2. In cither contradiction 
results since Ec{M) < by Definition 2.7 and U = {u\,u\Y satisfies the first- 
order system (2.2a). Hence, '&{z,x) is well-defined on [a;o,c]. For x € [xo,c) and 
M e 'D{z,c,xo,Oi), (3.2) implies that 

2ia{c, Xo, z){ui{z, x)*U2{z, x) — U2{z,x)*ui{z,x)) < 0. (3-11) 



(3.8) 
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This is equivalent to 

[ui{z,x)* - ia{c, xo, z)u2{z, x)*][ui{z, x) + ia{c, xq, z)u2{z, x)] 

< [ui{z,x)* + ia{c, Xo, z)u2{z, x)*][ui{z, x) - ia{c,xo, z)u2{z,x)] ^'^■^■^^ 

on [a;o, c). Given the nonsingularity of ui{z, x) — ia{c,Xo, z)u2{z,x) on [a;o, c], (3.12) 
imphes (3.9), with nonnegativity holding at x = c. 

Next, let M e C™^'", and suppose that '!?(z, x), defined by (3.8), is well-defined on 
[a;o,c], and satisfies (3.9). Then, on [xq,c), (3.12) and consequently (3.11) follow, 
which implies that (3.2) holds, and hence that M £ I?(2;, c, a;o, a). The analogous 
characterization of 'D±{z, xq, a) follows from Definition 2.12. □ 

By Lemma 3.1 one notes, for z e C\]R, that M e T^^z, c, xq, a) if and only if 

V{z,x,xo,a) = U2{z,x,xo,a)ui{z,x,xo,a)~^, xG[xo,c), (3.13) 
is well-defined while satisfying 

a{c,xo,z)lm{V{z,x,xo,a))>0, xG[xo,c). (3-14) 
In terms of V{z, x, xq, a) and by (3.8), one notes that 
'd{z,x,xo,a) = [Im + ia{c,xo, z)V{z,x,xo,a)] x 

(3.15) 

X [Im - icr{c, Xo, z)V{z, x, xo, a)] , x & [xo, c), 

is a Caylcy-type transformation of V"(z, x, xq, a). In the scalar context, this trans- 
formation corresponds to a conformal mapping of the complex upper half-plane to 
the unit disk. Moreover, defined as it is, V{z, x, xq, a) satisfies a Riccati differential 
equation that is associated with the Hamiltonian system (2.2a) while 'd{z,x,xo,a) 
satisfies a Riccati equation obtained by the Cayley-type transformation (3.15) ap- 
plied to the differential equation satisfied by V{z,x, xo, a). 

For the Dirac-type case of (2.2a), one observes by a simple calculation that 
V{z,x,xo,ao) is seen to satisfy a particular initial value problem for a Riccati 
differential equation. 

Lemma 3.3. Assume Hypotheses 2.1, 2.2, and 2.17. Let ao = (/„i 0) € C'"''^™^ 
let Uj{z,x) = Uj{z,x,Xo,cx.o), j = 1,2, be defined by (2.17), and suppose that 
V{z,x,xo,Oio) is well-defined by (3.13). Then, V{z,-) = V{z,-,xo,(Xo) satisfies, 

V'iz, x) + zV{z, xf + V{z, x)B2,2{x)V{z, x) + Bi^2{x)V{z, x) + V{z, x)B2,i{x) 
+ Bi^-l{x) + zim = 0, (3.16a) 
V{z,xo) = M, (3.16b) 

where Bj,k e C"^™, j,k = 1,2, are defined in (2.1d). 

Hence, by Lemma 3.1, the associated relations (3.13) and (3.14), and the unique- 
ness of solutions for (3.16), we obtain the following result for the Dirac-type case. 

Theorem 3.4. Assun^e Hypotheses 2.1, 2.2, and 2.17, and let qq = {Im 0) € 
(j^mx2m_ j'fi^fi^ M e I?(z, c, Xq, ao) if and only if the initial value problem given by 
(3.16) has a solution, V{z,-), well-defined and satisfying 

a{c, Xo, z)\m{V{z, x)) > 0, xe[xo, c). (3-17) 

Moreover, M e 'D±{z, xo, cto) if and only if (3.17) holds for c = ±oo. 
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Remark 3.5. An important consequence of Theorem 3.4 and the uniqueness of so- 
lutions for (3.16) is that solution trajectories for (3.16), which satisfy (3.17), consist 
of elements of Weyl disks; that is, 

V{z,x,XQ,ao) G 'D{z,c,x,aa), xG[xq,c). (3.18) 

Given the characterization of X'(z, c, xq, ao) in Defintion 2.7A, for each x G [xo,c) 
there is a /? e C^^^m ^j^j^ ^^^^ z)Im{P2pi) > 0, such that 

V{z, X, xo, ao) = M{z, c, x, ao, P). (3.19) 

It is in this sense that we let M{z, c, x, ao) denote our solution of the initial value 
problem (3.16) that satisfies (3.17). Analogously, 

V{z, X, xq, ao) e 'D±{z, X, ao), x € [xq, ±oo), (3.20) 

for trajectories of (3.16) that satisfy (3.17) for c = ±oo. Hence, in this sense, we 
let M±{z,x,ao) denote those solutions of (3.16) that satisfy (3.17) for c = ±oo. 
However, by Lemma 2.15, our Dirac system is in the limit point case at ±oo. 
Each 'D±{z, X, ao) consists of a unique matrix, and thus M±{z, x, ao) describes 
unique trajectories for (3.16a). This contrasts with the matrix-valued Schrodinger 
case considered in [20] where there are as many trajectories, each denoted by 
either M-^.{z,x,ao) or M-{z,x,ao), as there are matrices in a given initial disk 
T>±{z, xo,ao). 

Now for the Dirac-type case (2.4) with ao = (/„ 0) e C"'^^™^ ^^^^ ^z,x) = 
'&{z,x,xo,ao) defined in (3.8) and (3.15), and with x € [xo,c), one concludes that 

d{z, x)[ui{z, x) — ia{c, xo,z)u2{z, x)] = ui{z, x) + ia{c, xq, z)u2{z, x), (3.21) 

and hence that 

Im + ^{z, x) = 2u\{z, x)[u\{z, x) — ia{c, xq, z)u2{z, x)]^^, (3.22a) 

Im — '&{z,x) = —2ia{c, Xo, z)u2iz, x)[ui{z, x) — ia{c, Xo, z)u2{z, x)]~^ . (3.22b) 
Differentiating (3.21) one obtains 

{Ui - iaU2) = {Im - '&){ZU2 + B2,lUi + B2,2U2) 

(3.23) 

+ ia{Im + '>9){-ZUi- Bi^lUi- Bi^2U2)- 

By (3.22) one concludes that '&{z, ■, xq, ao) satisfies the initial value problem given 
by 



-ia{c,xo,z){zIm + Bi^i{x)) Bi^2{x) 

B2,i{x) ia{c,xo,z){zIm + B2,2{x)) 



X 



'd{z, Xo) = {Im + icric, xq, z)M){Im - i(7{c, Xo, z)M)-\ (3.24b) 

where Bj^k e C"^™, j. A: = 1, 2, satisfy Hypothesis 2.1. 

By Lemma 3.2 and the uniqueness of solutions for (3.24), one obtains the fol- 
lowing result in the Dirac-type case (2.4). 
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Theorem 3.6. Assume Hypothesis 2.17. Then M e 'D{z,c,xo,ao) if and only if 
the initial value problem given by (3.24) has a solution, '&{z, •) which is well-defined 
on [xo,c] and satisfies 

Im-Hz,xy§{z,x) >0, x€[xa,c). (3.25) 

Moreover, M e 'D±{z, xq, ao) if and only if (3.25) holds for c — ±00. 

Given the positivity present in (3.25), we note the exact correspondence which 
exists, by (3.15), between solutions of (3.16) that satisfy (3.17) and those sohitions 
of (3.24) that satisfy (3.25). In particular, given Remark 3.5, we rewrite (3.15) as 

'd{z,x,Xo,ao) = [Lm + i(T{c,XQ, z)M{z,c,x,aa)] x 

(3.26) 

X [Ijn - ia{c,xo,z)M{z,c,x,ao)] , xe[xo,c), 

Moreover, our Dirac system is in the limit point case at ±00. Consequently, there 
are unique solutions of (3.24), 'd±{z,-,xo,ao), z G C\IR, which satisfy (3.25) for 
c = ±00, and which correspond to the unique solutions of (3.16), M±{z,x,ao), 
which satisfy (3.17) for c = ±00; specifically, 

'3±{z,x,Xo,ao) = [Im ±i(^{z)M±{z,x,ao)][ImT i(^{z)M±{z,x,ao)]~^ ■ (3.27) 

These relationships form the basis for the analysis to follow. The asymptotic 

result (3.1) is obtained by an analysis of the corresponding asymptotic behavior 
for all solutions '&{z,-,Xo,ao) described in (3.24), these include among them the 
particular solutions '&±{z,-,xo,ao). Thus asymptotic behavior is deduced for all 
corresponding solutions M (z, c, •, ao) of (3.16) which include among them the solu- 
tions M±{z, •, ao)- The advantage of this approach comes from the compactification 
inherent in the Cayley-type transformation (3.26), and the resulting boundedness 
of the solutions as a consequence of (3.25). 

We pause for a moment to address, in the following remark, a point raised by us 
in [20] for the matrix- valued Schrodinger case described in (2.5). 

Remark 3.7. With Uj{z,x) = Uj{z,x,xo,C(), j = 1,2, defined in (2.17) for the 

general Hamiltonian system (2.2a), an analog to Lemma 3.2 for the characterization 
of "D^z, c, xo, a) is obtained by replacing the expression in (3.7) with 

ui{z, x) - i\z\~^^'^a{c, Xo,z)u2{z, x), (3.28) 

and by replacing the definition for i!}{z, x) = •&{z, x, Xq, a) given in (3.8) with 

t?(z, x) ={ui{z, x) + i\z\~^l'^a{c, xo, z)u2(z, x)) x ^ 
X (u\{z,x) - i\z\~^l'^a{c,XQ,z)u2{z,x)Y^ . 

Specific to the matrix-valued Schrodinger case, we obtain analogs of Lemma 3.3, 
Theorem 3.4, and Theorem 3.6 by replacing equation (3.16a) with 

V'[z, x) + V{z, xf - Q{x) + zlm = ^ (3.30) 

and by replacing the equations in (3.24) with 

.',,.)=„(c...,4(;-:|-);)'(-N-(^;™-«w) ^,^|_o,^_^ 

■&{z,xq) = {Im+i\z\-'^^'^(T{c,xo,z)M){Im-i\z\-'^/'^a{c,xo,z)M)-\ (3.31b) 
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^'{z,t) = \\z\-' 



(3.34a) 



V^iz, c, xo, ao) was defined in [20] to be the set of those M e C™^™ for which the 
intial vahic problem given by (3.31) has a solution, ??(z, x), which is well-defined on 
[a;o, c] and satisfies (3.25). In [20] we showed that V^z, c, xq, ao) C 'D^{z, c, xo, cxo)- 
This was sufficient for the subsequent analysis in [20]. However, as the analog of 
Theorem 3.6 now shows, one actually has equality of the two disks in [20], that is, 

'D{z, c, Xo, ao) = V^{z, c, xq, ao). (3.32) 

To obtain a proof of (3.1) for the Dirac-typc case, wc adapt an approach due 
to Atkinson [12] for proving a result analogous to (3.1) for the matrix- valued 
Schrodinger case (cf., e.g., [20, Theorem 3.1]) In light of Remark 3.12, we begin 
by restricting our attention to ^ e C+, and as in the previous discussion, take 
ao = {Im 0) e C^x^'". 

First we introduce two systems related to (3.24) by means of a change of vari- 
ables. Let 

ip{z,f) = ■d{z,x), t={x — xo)\z\, x€l[xo,c). (3.33) 
With this change, (3.24) becomes 

l,^,_i /^/„ + ^(z,t)*y /T«(^/m + Si,i(i)) Bi^t) \ 

With M = M{z, c, Xo, chq) G 'D{z, c, xo, chq) (3.24b) becomes 

ip{z, 0) = {ilm T M{z, c, Xo, ao)){iIm ± M{z, c, xq, ao))"S (3.34b) 
and (3.25) becomes 

<fi{z, t)*<p{z, t)<Im t€ [0, (c - xo)\z\), (3.34c) 

where in (3.34a), 

BjMt) ^ Bj^kixo+t\z\-^), j,k = 1,2. (3.34d) 

In the complete system (3.34), one now has a set of conditions equivalent to system 
(3.24) and (3.25). 

We recall that Cg C C+ represents the open sector with vertex at zero, symmetry 
axis along the positive imaginary axis, and opening angle e, with < £ < 7r/2. Next, 
consider a sequence, 2;„ e C^, n e N, such that |z„| ^ oo as n — > oo and such that 

< e < 5n = arg (zn) <7v — s. (3.35) 

By choosing an appropriate subsequence, we may assume that 

6n^S €[5,^-8]. (3.36) 

Let (f{zn,t) denote a corresponding sequence of functions that satisfy (3.34a) and 
(3.34c), with initial data, ^p{zn,Q), defined by (3.34b) for a sequence of points 
M{zn, c, Xo, ao), where each M{zn, c, xq, uq) is chosen to be an element of the disk 
'D{zn,c,XQ,ao). Note that as z„ oo, the intervals described in (3.34c) even- 
tually cover all compact subintervals of IR.+ . Given the uniform boundedness of 
'■Pn{t) = ^{zn,i) described in (3.34c), we assume, upon passing to an appropriate 
subsequence still denoted by iy5n(0), that 

<Pn(0) = <p(2;„,0) V±{'5), for ± (c— xo) > as n ^ oo, (3.37) 
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and as a consequence, that 

^±{S)*<p±{5) < Im. (3.38) 

With (p±{S) defined in (3.37) as |2;„| ^ oo, we consider limiting systems associ- 
ated with (3.34): 

(3.39a) 

V±{0) = ^±iS). (3.39b) 

Theorem 3.8. Assume Hypothesis 2.17. Then the solution ri± of (3.39) satisfies 

V±{trv±{t) <Im, tG [0, ±oo). (3.40) 

Moreover, the solutions ipn = ^{zn, ■) of (3.34) converge to r}± uniformly on [0, ±T] 
for every T > 0, as n — > oo. 

Proof. In this proof, wc consider only the case corresponding to t > 0, that is, 
7y+(0) = tp+{S) in (3.39b). The other case foUows in a similar manner. For this 
reason, we let r]{t) = r]+{t) in the remaining discussion. We also let T e M+ be the 
greatest value such that (3.40) holds for t G [0,T] and show that (3.40) must hold 
for some [0, T'] with T' > T, thus proving T = oo. 

The solution of (3.39), r], presumed to be defined on [0, T], can be continued onto 
some [0, T'] with T' > T; r] then satisfies 

V{t)*v{t) < t^^Im (3.41) 

for <t <T' and for some k > 1. 

For brevity, let iy?J^(t) = Gn{fn,t) denote (3.34a) with z = z„, and let r]'{t) = H{r], t) 
denote (3.39a) in the following. Integrating (3.39a) and (3.34a), one obtains 

Vn{t) - v{t) = Vn{0) - foiS) + / ds{Gn{v, s) - H{r], s)} 

Jo 

+ [ ds{Gn{'Pn,s)-Gn{v,s)}. (3.42) 

Jo 

We note that 

Gn{v,s) - H{r^,s) = \i{e'' - e^^")(/„ + r,(s))2 - iz(e^^ - e^^")(/„ - in{s)f + 

2 

+ Y,F3,k{zn,s), (3.43) 



where. 



Fi,i{zn,s) = ~i\zn\-^{Im+'n{s))B-i,i{s){Im+ri{s)), (3.44a) 

F2,2{Zn,s) = ^\Zn\-\lm-r]{s))B2,2{s){Im-r]{s)), (3.44b) 

Fl,2{Zn,s) = ]^\Zn\-\l^ + v{s))BiaWrn - riis)), (3.44c) 

i=^2,i(^„, s) = ^\zn\-\lm ' r?(s))B2,i (s)(/„, + r?(s)). (3.44d) 
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Thus, for t e [0, T'], (3.41) implies that as n ^ oo 



10 

and together with (3.33) and (3.34d) that 

ft / /.£i;o+f|z„r^ 



/ \\Im±v{sWc^.r.ds = o{l), (3.45) 
Jo 



yjl^i.fc(s)llc-x^rfs = 0(^y ||B,-fc(s)||c™x™rfsj =0(1). (3.46) 

(Here || • Hcnx™ denotes a norm on C™^™.) Hence, by (3.43)-(3.46), one infers that 
for t e [0, T'] and as n — > oo, 

rt 



[ {Gn{v,s)-H{r,,s)}ds = o{l). (3.47) 
Jo 



Next, one notes that 



2 

Gn{<Pn,s) - Gn{v,s) = 2ie'^- {r]{s) - <fin{s)) + Kj^k{Zn,s), (3.48) 

j,k=l 

where 

Kl,l{Zn,s) = —\Zn\~'^{{Im + 'Pn)Bl,l{s){ifin - V) + {^n - V)Bl,l{s){Im + V)} , 



K2,2{Zn,s) = -\Zn\ ^ {{Im - fn)B2,2{s){ll - fn) + {V - fn)B2,2{s){Im - V)} 
Ki,2{Zn,s) = ^l^^nPHUm + <Pn)-Bl,2(s)(r? - fn) + {^n - '7)-Bi,2 (s) (/m " V)} 



K2,l{Zn,s) = l-\Zn\ ^ {{Im - fn)B2,l{s){(pn - V) + iV - <fn)B2,l{s){Im + V)}- 



(3.49a) 

?)}, 

(3.49b) 

?)}. 

(3.49c) 

?)}• 
(3.49d) 

By (3.38) and (3.41), for s e [0,T'], 

\\Im±Vn{s)\\c"^xm <2, ||7„±?7(s)||cmxm < K+1, (3.50) 

and hence by (3.48)-(3.50), 

\\Gn{iPn, S) - Gniv, s)||c"'X'" 

< Us)-^^{s)\\c^.^{2+\^{3 + k) J2 WBMWc^^-]- (3.51) 

Of course, by (3.37) as n — > oo, 

<A„(0)-<^+(^)=o(l). (3.52) 
Thus, by (3.46), (3.51) and (3.52), one concludes for t e [0, T'] and as n ^ oo, that 
lkn(i)-r?(i)||c— <o(l) 

+ £ll^n(s)-r?(s)||c™x™|2+^I^(3 + «) ^^ |l4fe(s)|lc™x.|ds. (3.53) 
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Gronwall's inequality applied to (3.53) together with a consideration of the effect 
of the variable change (3.33), as illustrated in (3.46), yields 

'Pnit) - ri{t) ^ as n ^ oo (3.54) 

uniformly for t S [0,T']. Thus by (3.34c), the contradiction results that for all 

t e [0, r'], 7] satisfies (3.40). □ 

What solutions of (3.39) satisfy (3.40)? 

Lemma 3.9. Assume Hypothesis 2.17. Ifrj± is a solution of (3.39a) which satisfies 
(3.40), then 

= 'n±{t), t e [0, ±oo). (3.55) 

Proof. We note that (3.39a) is equivalent to (3.34a) with B = Q. By the vari- 
able change (3.33), (3.39a) is also equivalent to (3.24a) with B = Q. Next, we 
recall the connection between the Riccati-type equations (3.24a), and (3.16a) by 
means of the Cayley transformation (3.26). Solution matrices of (3.39a) which 
statisfy (3.40) &i t — Q thus correspond to solution matrices, V{z,-), of (3.16a) 
for which lui(y{z,XQ)) > 0. Moreover, solutions of (3.16a) for which for which 
lm{V{z,xo)) > are obtainable from solutions of (2.2a), with B = 0, by means of 
(2.17) with Im(i\/) > 0. Thus, by utilizing this connection between explicit expo- 
nential solutions of (2.2a) with B = and solutions of the Riccati-type equation 
(3.16a), and by performing on the resulting solution of (3.16a) the conformal map- 
ping (3.26) followed by the variable transformation (3.33), one obtains the following 
solution for (3.34a), 

^iz,t) ^ {iI^TM){iI^±M)-^cxp{T2ite'^), (3.56) 

for ±t > 0, Im(±M) > 0, and z G C+. By hypothesis, < 6 < n. Thus the 
exponential term in (3.56) will result in 

\\ip{z,t)\\cmxm > 1 ast^ioo (3.57) 

unless 

M = ±ilm, (3.58) 
thus implying (3.55). □ 
One then obtains the following result. 

Corollary 3.10. With (t>±{S) defined in (3.37), rj±{0) = (j)±{6) = 0. 

For M{zn,c,xo,ao) G 'D{zn,c, xo,(Xo), it follows by (3.34b), (3.37), and Corol- 
lary 3.10 that 

[HmT M{zn,c,xo,ao)][iIm± M{zn,c,xo,ao)]~^ =o(l), ±{c-xo) > 0, 

(3.59) 

as n ^ oo. Hence one infers, for elements of ^{zn, c, xq, ao), that 

M{z„,c,xo,ao) =±iIm + o{l), ±(c-a;o)>0, (3.60) 

as \z\ ^ oo in C^. This proves (3.1). Actually, (3.60) is a statement for all elements 

of 2^(2, c, xq, ao) including the particular element AI±{z, xq, ao), for ±(c — xq) > 0. 

In (3.1) an asymptotic expansion is given that is uniform with respect to arg(2;) 
for 1^1 00 in Ce. We now vary the reference point, xq, and observe that the 
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asymptotic expansion in (3.1) is also uniform with respect to xq whenever xo is 
confined to a compact subset of M. 

Theorem 3.11. Assume Hypothesis 2.17. Let ao = {Im 0) € C™''^™^ denote 
by Ce C C+ the open sector with vertex at zero, symmetry axis along the positve 
imaginary axis and opening angle s, with < s < ■jt/2. Let M±{z,xo,ao) be the 
unique elements of the limit disks 'D±{z, xq, ao) for the Dirac system given by (2.2) 
and (2.4). Then, 

M±{z,x,ao) = ±iIm + o{l) (3.61) 

zeCe 

uniformly with respect to arg(^;), for \z\ ^ oo in C^, and uniformly with respect to 
x, as long as x varies in compact subsets of [a;o;±oo). 

Proof. We note that the system (3.39) is independent of the reference point xq. 
Next, we recall that 6, defined in (3.36) is determined by an apriori choice of 
the sequence z„, subject only to z„ being in Cg (c.f. (3.35)). Moreover, wc note 
that <f±{6), defined as a limit in (3.37), described explicity in Corollary 3.10, and 
which gives solutions of (3.39) satisfying (3.40) for t G [0, ±oo), is also independent 
of the reference point xq. Thus, had wc chosen a different point of reference, 
Xq ^ Xq, at the start, the asymptotic analysis begun in Theorem 3.8 and continued 
through (3.59), would remain the same after the variable change in (3.33), except 
for the integral expression in (3.46) in which xq would be replaced by Xq. However, 
given the local integrability assumption on B in Hypothesis 2.1, one concludes that 
this integral expression is uniformly continuous with respect to xq whenever xq 
is confined to a compact subset of M. Thus (3.46), and consequently (3.54), are 
uniform with respect to t and with respect to xq whenever both are confined to 
compact subsets of M. Consequently, (3.59) holds for elements T>{z,c,xo,ao), that 
this asymptotic expansion is uniform with respect to arg(z) for |z| ^ oo in C^, and 
that it is uniform with respect to xq when xq is confined to compact subsets of 
M. □ 

Remark 3.12. (i) In the special case m = 1, the leading-order asymptotics (3.61) 
was published by Evcritt, Hinton, and Shaw [32] in 1983. For asymptotic estimates 
of Weyl solutions in the case m = 1 we refer to [97]. 

(ii) A comparison of (3.61) with (2.60) then proves that the leading-order asymp- 
totic behavior (3.61) is in fact independent of the boundary condition at xq indexed 
by a, that is, 

M±(z,a;o,a) = ±i/„ + o(l) (3.62) 

zECe 

for any a satisfying the conditions stated in (2.9). In the scalar case m = 1 this fact 
had been noticed in [32]. This boundary condition independence of the leading- 
order asymptotic behavior of M± {z,Xo,a) is in sharp contrast to the case of matrix- 
valued Schrodinger operators (see, e.g., [20]). Moreover, regarding the conclusion 
of Theorem 3.11, no generality is lost by assuming that C C-|- because of (2.57). 

4. Higher Order Terms in the Asymptotic Expansion of M±{z,x,a) 

In this section wc shall prove one of our principal results of this paper, the as- 
ymptotic high-energy expansion of M^{z, x, ao) to arbitrarily high orders in sectors 
of the type C C+ as defined in Theorem 3.11. 
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Throughout this section we choose z e C-|- . We also recall the following notion: 

X € [a, 6) (rasp., x G (a, 5]) is called a right (resp., left) Lcbesgue point of an element 
q G L^{{a, b)), a < b ifjj dx' \q{x+x') — q{x)\ = o{e) (resp., /J dx' \q{x~x') — q{x)\ = 
o(e)) as e J, 0. Similarly, x G (a, 6) is called a Lebesgue point of g G L^{{a,b)) if 
He ^ '?(■'')! = o(£) as £ J, 0. The set of all such points is then denoted 

the right (resp., left) Lebesgue set of q on [a, b] in the former case and simply the 
Lebesgue set of q on [a,b] in the latter case. The analogous notions are applied 
to 2m X 2m matrices B G L^((a, 6))^'"^^'" by simultaneously considering all 4to^ 
entries of B. The right (resp., left) Lebesgue set of B on [a,b] is then simply the 
intersection of the right (resp., loft) Lebesgue sets of Bj^k for all 1 < j, k < 2m, and 
similarly for the Lebesgue set of B, etc. 

Finally, wc need one more ingredient, recently proven by Rybkin [101, Lemma 3] 
using appropriate maximal functions. Let q G L^{{xo,oo)), supp((/) C [xo,Xo + R] 
for some R> 0, and suppose x G [xo,xo + R] is a right Lebesgue point of q. Then 

dx' g(x') exp(2iz(x' - x)) = -^ + o{\z\-^). (4.1) 



L 



An alternative proof of (4.1) follows from [117, Theorem 1.13], which implies 

lim dx'\q{x')-q{x)\eyiY>{2iz{x' -x))=Q (4.2) 

kHoo Jx 

for any right Lebesgue point x of q. 

We start with the simpler case where B has compact support contained in some 
interval [xQ,yo]. Below in (4.3) and in analogous formulas in this section, || • ||c<x« 
denotes a norm in 

Lemma 4.1. Fix Xq, yo G M with yo > xq and let x > xq- Suppose A = l2m, B G 
L^{[xa,xo -I- i?])2'nx2m all R > 0, B = B* a.e. on (xq.oo). In addition, assume 
that B has compact support contained in [xo,yo], that G {[xq, yo])'^"^^'^"^ 

for some G N, that x is a right Lebesgue point of B^^~^\ and that 



ess sup 

ye[xQ,vo] 



|z|— »oo 

zeCe 



/ dx' B^^-'\x')eM^iz{x' -y)) + —B^^-^\y) 
o{\z\-'). (4.3) 



If N = 1, suppose in addition Bk^k'Bgj, G L\[xo,yo])"'''"' for all k,k',i,i' G 
{1,2}. Let ao = {Im 0) G C^^^™ and denote by M+{z,x,ao), x > Xq, the unique 
Weyl-Titchmarsh matrix associated with the half-line Dirac-type operator D+{ao) 
in (2.88). Then, as \z\ — > oo m C^, M+{z,x,ao) has an asymptotic expansion of 
the form 



N 



M+{z,x,ao) = iIrrt + y"^m+^k{x,ao)z-'' + o(\z\-^), N eN. (4.4) 



zeCe 



fe=l 



The expansion (4.4) is uniform with respect to arg (z) for \z\ — > oo in and 
uniform in x as long as x varies in compact subintervals of [xq, oo) intersected with 
the right Lebesgue set of B^^~^h The expansion coefficients m+,fe(a;, ao) can be 
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recursively computed from 

1 i 
m+^i{x,ao) = --{Bi^2{x) + -82,1(0;)) + -{Bi^i{x) - -62,2(2;)), 

■ / k 
m+,fe+i(a;, ao) = - I m'_^_^f,{x, ao) + ^ m+,^(a;, ao)m+,fe+i_^(a;, ao) 
^ 1=1 

k 

+ ^ m+^e{x, ao)B2,2(a;)m+,fe_^(a;, ao) (4.5) 

+ Bi^2ix)m+^k{x, ao) + m+,fc(a;, ao)B2,i{x) 

l<k<N-l. 

Proof. In the following let z G C_|_, and x > xq. The existence of an expansion of 
the type (4.4) is shown as follows. First one considers a matrix Volterra integral 
equation of the type 

U+{z,x,aQ) = ^l^p^ exp{iz{x — Xq)) + J dx'K{z,x,x')JB{x')U^{z,x',ao), 

(4.6) 

where 

[?4.,.,«o)^(?-;<;::^;:j)ci'(|...-))--, (4.7) 

and K abbreviates the 2m x 2m. Volterra Green's kernel 

, M _ / cos(2;(x - x'))/m sm{zix- x'))Im\ s^ 
i^[z,x,x)- ^_ ^.^^^^^ _ ^^^^^^^ _ . ^4.»; 

Clearly, U+{z,-,aQ) solves the Dirac-typc system (2.2) and (2.4). In addition, it 
satisfies U+{z,-,ao) € I/^([a;o, oo))^"*^^"*. Thus, up to normalization, C/+(z,-,ao) 
represents the Weyl solution associated with B on the half- line [a;o,oo). Next, 
introducing 

V^{z,x,ao) = =U+{z,x,ao)eM-i4x-xo)), (4.9) 



^v+^2{z,x,ao) 
one rewrites (4.6) in the form 



V+{z,x,ao) = [J^^ + J\x' K{z,x,x')JB{x')V+{z,x',ao), (4.10) 
K{z, X, x') 



where 

1 / (1 + exp(2iz(a::' — a;)))/m — i(l — exp(2iz(x' — a:)))/m\ 

2 li(l — exp(2i2;(a;' — x)))Im (1 + exp(2i^(a;' — x)))Im I 

(4.11) 

Thus, one infers, 

M+{z,x,ao) = u+^2{z,x,ao)u+^i{z,x,ao)~^ = v+,2{z,x,ao)v+^i{z,x,ao)~'^ . 

(4.12) 

Introducing 





-tC2\ 




lD2\ 


[iCi 


C2 j ' 




D2J 



R= A/ I S= - , (4.13) 
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where 



<^l = -Sl*,2-«Sl,l. ^2 = 5i,2 - 1^2,2, (4.14) 

Di = -BU + i-Bi,i, £»2 = -Bi,2 + iB2,2, (4.15) 



(4.10) results in 

V+{z, X, ao) = (^Ip ^ + dx' {R{x') + S{x') exp{2iz{x' - x)))V+{z, x', ao) 



(4.16) 



2"'' / dxi {R{xi) + 5(a;i)e2^^(^i-^)) x 
fe=i "'^^ 

X / dx2(i?(x2) + 5(x2)e2*"(^^-"i))... (4.17) 
rVo 

■■■ dxk {R{xk) + S{xk)e''"^="'-'"'-' 



il„ 



This yields 

\\v+j{z,x,ao)\\<Cj, z€C+,lm{z)>0,x>xo,j = l,2 (4.18) 

for some Cj > 0, j = 1,2, depending on Integrating by parts in (4.17), 

repeatedly applying (4.1) and (4.3) to q{x) = {S{x))j^k for all 1 < j, A: < 2m then 
results in the existence of an asymptotic expansion for V+(z, x, ao) of the type 



y+(.,.,ao)= =E^'+.^(-'"o).-'= + o(N-")^ (4.19) 



Inserting the expansions for v^^2{z,x,aQ) and v^^i(z,x,a())^^ into (4.12) (using 
a geometric series expansion for v+^i{z, x,ao)~^) then yields the existence of an 
expansion of the type (4.4) for M^{z, x, ao). The actual expansion coefficients and 
the associated recursion relation (4.5) then follow upon inserting expansion (4.4) 
into the Riccati-type equation (3.16a). The stated uniformity assertions concerning 
the asymptotic expansion (4.4) then follow from iterating the system of Volterra 
integral integral equations (4.10). □ 

Remark 4.2. The analogous solution U-{z, -jao) of the Dirac-type operator (2.75) 
on the interval (— oo,a;o] satisfies 

U-{z, x,ao) = ) exp{-iz{x - xo)) 



J — c 



dx' K{z, X, x')JB{x')U- {z, x', ao), (4.20) 



with integral kernel K given by (4.8). (Again U- coincides with the Wcyl solution 
U- up to normalization.) A closer look at the system of Volterra integral equations 
(4.6), (4.16), (4.17), and similarly in connection with (4.20), then reveals that 
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U±{z,-,ao) have the asymptotic behavior 



N 



(4.21) 

with leading asymptotics determined as follows. 

v±,Q,i{x,ao) = Im + w±fi^i{x,ao), 

~ ~ (4.22) 

^^±,o,2(a;,ao) = ±i{lm +w±fi^i{x,aQ)), 

where w5±,o,i(a;, ao) satisfies 

'u;±^o,i(^'"o) = ^[-62,1(2;) - -61,2(2;) ±iB2,2ix) ±iBi^i{x)] (/^ + w±,o,i{x,ao)), 

(4.23) 

and 

lim w± i{x, ao) = (4.24) 

(in fact, v±fi^i{-,ao) = Im, v±,o,2{j'^o) = ±ilm, and v±^k,j{-,ao) = 0, j = 1,2, 
1 < k < N outside the support of B). In particular, 

u;±,o,i(a;,ao) = (4.25) 

and hence 

U±{z,x,ao) , = ((J"! ) +0(1)) exp(±zz(x-xo)), (4.26) 

zeCe 

if and only if B is in the normal form 

B{x) = [ f 'ij^j J ' 'Br,i(x) = B,,,ix), Bl^ix) = B,,2{x) a.e. 

\^l,2(X) —tSi^i{X)J 

(4.27) 

For more details we refer to Lemma 5.1. 

Next we recall an elementary result on finite-dimensional evolution equations 
essentially taken from [98] (cf. also [20, Lemma 4.2]). 

Lemma 4.3. ([98].) Let Tj e Lj^^(M)'"^™, j = 1,2. Then any m x m matrix- 
valued solution X of 

X'{x) = Ti{x)X{x) + X{x)T2{x) for a.e. a; e M, (4.28) 

is of the type 

X{x) = Y{x)CZ{x), (4.29) 
where C is a constant mxm matrix and Y is a fundamental system of solutions of 

^'{x)=Ti{x)^{x) (4.30) 
and Z is a fundamental system of solutions of 

^'{x) = ^x)T2{x). (4.31) 
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The next result provides the proper extension of Lemma 4.3 in [20] in the context 

of matrix- vahicd Schrodinger operators (which in turn extended Proposition 2.1 in 
the scalar context in [49] to the matrix-valued case) to the Dirac-type case under 
consideration. 

Lemma 4.4. Fix a;o,yo S K with yo > xq. Suppose Aj = l2m, Bj e L^{[xo,xo + 

^])2mx2m ^ > Q, Bj = B* a.e. on [xo,oo), j = 1,2, and Bi = B2 a.e. on 

[xo,yo]- Let ao = {Im 0) € C"^^"* and denote by Mj^+{z,x,ao), x > xo, the unique 
Weyl-Titchmarsh matrix corresponding to the half-line Dirac operators D+j{ao), 
j = 1,2, in (2.88). Then, 

[M[^^{z,x,ao) - M2 _^(z,x,q;o)] 

= -{z/2)[Mi^+{z, X, ao) + M2,+ {z, x, ao)] [Mi,+(z, x, ao) - M2,+ {z, x, ao)] 

- {z/2)[Mi^+{z,x,ao) - M2,+{z,x,ao)][Mi^+{z,x,ao) + M2^+{z, x, ao)] 

- [Mi^+{z,x,ao) + M2,+ iz,x,ao)]B2,2ix)[Mi^+{z,x,ao) - M2,+ (z, x, q;o)]/2 

- [Mi^+{z,x,ao) - M2,+ {z,x,ao)]B2,2{x)[Mi^+{z,x,ao) + M2,+ {z,x,ao)]/2 

- Bi^2{x)[Mi^+{z, X, ao) - M2,+ {z,x,ao)] 

- [Mi^+{z,x,ao) - M2^+{z,x,ao)]B2^i{x) for a.e. x G [xo,yo], (4.32) 
where we denoted Bi = B2 = (^bI'I B22) (^o,2/o)- 

Proof. This is obvious from (3.16a). □ 

Lemma 4.5. Fix xo,yo G ^ with yo > xo- Suppose Aj = l2m, Bj G L^{[xo,xo + 
R])2mx2m ^ > 0, and Bj = B*j a.e. on [xq, 00), j = 1, 2. Let ao = {Lm 0) S 

Qmx2m denote by Mj_+(z,x,ao), x > xo, the unique Weyl-Titchm,arsh w,atrix 
corresponding to the half-line Dirac operators D^_j{ao), j = 1,2, in (2.88). Define 

Ti{z,x) = -{zl2)[Mi^+{z,x,ao) + M2,+ {z,x,ao)] 

- (l/2)[Mi,+ (2, X, ao) + M2,+{z, x, ao)]B2,2(x) - Bi,2(a;), (4.33) 
T2{z, x) = -{z/2)[Mi^+{z, X, ao) + M2,+{z, x, ao)] 

- {l/2)B2,2{x)[Mi,+{z,x,ao) + M2,+ {z,x,ao)]-B2,i{x), (4.34) 

for a.e. x € [xo,yo]- In addition, assume Y+{z,-) and Z+{z,-) to be fundamental 

matrix solutions of 

'^'{z,x) =Ti{z,x)'^{z,x) and <^'{z,x) = '^{z,x)T2{z,x) (4.35) 
on [xo,yo], respectively, with 

yo) = Im, Z+{z, yo) = Im- (4.36) 

Then, as \z\ 00, z G C^, 

||F+(2;,a;o)||c"'X"',||-^-K(-2,a;o)||c-x- < exp(-Im(2;)(2/o - a;o)(l + o(l))). (4.37) 
Proof. Define rj{z,x), j = 1,2, by 

{z,x)= Tj {z, x) + izim, i = 1, 2, (4.38) 

then 

/ dx]\rj{z,x)]\c^.^ = o{z), j = l,2 (4.39) 
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due to the uniform nature of the asymptotic expansion (3.61) for x varying in 
compact intervals. Next, introduce 

E+{z, X, yo) = Im exp{iz{yo - x)), x < yo, (4.40) 

then 

Y+{z,x)=E+{z,x,yo)- dx' E+{z,x,x')Ti{z,x')Y+{z,x'), (4.41) 

J X 

Z+{z,x) = E+{z,x,yo) - dx' Z+{z,x')T2{z,x')E+{z,x,x'). (4.42) 

J X 

Using 

||£;+(z,a;o,2/o)||c-xm < exp(-Im(2;)(i/o - ^o)), (4.43) 
a standard Volterra-type iteration argument in (4.41), (4.42) then yields 

||y+(2,a;o)||c-x„ <exp(^-lm(z)(yo-Xo) + dx\\fi{z,x)\\^ , (4.44) 

||Z+(^,a;o)||c-x„. <exp(^-Im(^)(yo-a;o) + ^'"'da;||r2(^,a;)||^ , (4.45) 

and hence (4.37). □ 

Theorem 4.6. Fix xq, j/o G R with yo > xo- Suppose Aj = l2m, Bj e L^{[xo, xo + 
R]frnx2m /J > Q, Bj = B* tt.e. ou [ajQ, oo), j = 1,2, and Bi = B2 a.e. on 

[xo,yo]- Let ao = {Im 0) G ([]"ix2ni Q,^^ denote hy Mj,-i^{z,x,ao), x > xq, the unique 
Weyl-Titchmarsh matrix corresponding to the half-line Dirac operators Dj^j{ao), 
j = 1,2, in (2.88). Then, as \z\ — > 00 in C^, 

||Mi_+(2;,a;o,ao) - M2,+(2;,a;o,Q!o)||c"'xm < Cexp(-2Im(^;)(yo - a;o)(l + o(l))) 

(4.46) 

for some constant C > 0. 

Proof. Define for z G C\M, x e [a;o,?/o]) 

X+{z,x) = Mi,+{z,x,ao) - M2,+ {z,x,ao), (4.47) 
and for z e C\M and a.e. x G [a;o,yo]j 

Ti{z,x) = -{z/2)[Mi^+{z,xo,ao) + M2^+{z,xo,ao)\ 

- (l/2)[Mi,+ (0, Xo, ao) + M2,+ {z, xo, ao)]B2,2{x) - ^1,2(2;), (4.48) 
T2iz,x) = -{z/2)[Mi^+{z,xo,ao) + M2,+ {z,xo,ao)] 

- {l/2)B2,2{x)[Mi,+{z,xo,ao) + M2,+{z,xo,ao)] - B2,i{x). (4.49) 
By Lemma 4.4, 

= TiX+ + x+r2 (4.50) 

and hence by Lemma 4.3, 

X+{z,x) = Y+{z,x)X+{z,xi)Z+{z,x), (4.51) 
where Y^{z, x) and Z^{z, x) arc fundamental solution matrices of 

^'{z, x) = ri(z, x)^{z, x) and ^'{z, x) = $(z, a;)r2(z, x), (4.52) 
respectively, with 

Y+{z, yo) = Im, Z+{z, yo) = Im- (4.53) 
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By Lemma 4.5, 

||y+(2,a;o)||c"'xm, \\Z+{z,XQ)\\cmxm < exp(-Im(^;)(yo - a;o))(l + o(l))) (4.54) 
as \z\ oo, z G Cg. Thus, as 1^1 — > oo, ^ e Ce, 

||X+(2;,a;o)||c'"X'" < H^+C-z, yo)||c"'X"' ||i + (2^,a;o)||c"'X"' ||-^+(2^,a;o)||c"»xm 

< Cexp(-2Im(2)(yo - a;o)(l + o(l))) (4.55) 
for some constant C > by (3.61), (4.51), and (4.54). □ 

Given these preparations we can now drop the compact support assumption on 
B in Lemma 4.1 and hence arrive at one of the principal results of this paper. 

Theorem 4.7. Fix xg, yo G IK. with uq > xq and suppose A — /2m; B G L^{[xq, xo + 
j^j-j2mx2m y^^, R> 0, and B = B* a.e. on {xo,oo). In addition, assume that for 
some N gN, B^^''^^ e Z/^([a;o, c])^"^^'" for all c > xo, that xq is a right Lebesgue 
point of B'^^~^\ and that 



ess sup 



fVo 1 

/ dx' (x') exp(2iz(x' - y)) + B^^-^^ (y) 

Jv 2i2: 



C2mx2r> 



= o{\z\-'). (4.56) 

\z\ — ^oo 

If N = 1, suppose in addition Bk,k'Bi^ir e L^([a;o, yo])™^™ for all k,k',£,£' G 
{1,2}. Let ao = (/,„ 0) G ([^rnx2m denote by M^{z,Xo,ao) the unique element 
of the limit disk 2?+(z, .tq, ao) for the half-line Dirac operator D_|_(ao) in (2.88). 
Then, as \z\ oo in C^, M^{z, xo,ao) has an asymptotic expansion of the form 

N 

M+{z,xo,ao) = iIm + y^m+,k{xo,ao)z-'' +o(\z\-^), N gN. (4.57) 



k=l 



The expansion (4.57) is uniform with respect to arg(z) for \z\ oo in C^- The 
expansion coefficients m+,fc(a;o, ao) can be recursively computed from (4.5). 



Proof. Define 



B{x) = l^^^^ for a; e [a;o,yo], a;o < yo 
1 otherwise 



and apply Theorem 4.6 with Bi ~ B, B2 = B. Then (in obvious notation) 

\\M+{z, Xo, ao) - M+{z, xo, ao)||c™x™ < Cexp(-2Im(2;)(2/o - a;o)(l + o(l))) 

(4.59) 

as \z\ 00, z G Ce, and hence the asymptotic expansion (4.4) for M+(z, a;o, ao) in 
Lemma 4.1 coincides with that of Mj^{z, a;o, ao). □ 

In analogy to Theorem 3.11, the asymptotic expansion (4.57) extends to one 

for M+(z,a;, ao) valid uniformly with respect to x as long as x varies in compact 
subintcrvals of [a;o,oo) intersected with the right Lebesgue set of B^^^^K 

Theorem 4.8. Fix Xo G R and let x > xo- Suppose A = l2m, B G L^([a;o,a;o + 
^j^2mx2m y^j, i? > 0, and B = B* a.e. on {xo, oo). In addition, assume that for 
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some N gN, B(-'^-i) e I,^([a;o, c))^™^^™ for all c > xq, that x is a right Lebesgue 
point of B'^^~^\ and that for all R> 0, 

fXo+R 



ess sup 

ve[xo,xo + R] 



/ dx' ^(^-1) (a;') exp(2iz(x' - y)) + — B^^-^^ (y) 

Jy 



(J^2Tnx2T? 



= o{\z\-'). (4.60) 

\z\ — *oc 



If N = 1, suppose in addition Bk.k'Bgj' S L'^{[xo,xo + i?])™^'" for all R > 
all k,k',e,e' e {1,2}. Let oa = [Im 0) G C"^^™ ^^^^^g M+{z.x,ao), 
X > Xq, the unique element of the limit disk 'D+{z, a;, ao) for the half-line Dirac op- 
erator D^{ao) in (2.88). Then, as \z\ oo in Cg, M+(z,a;, ao) has an asymptotic 
expansion of the form 

N 

M+{z,x,aQ) = ilm + y" m+^k{x,ao)z-^ + o{\z\-^), TV e N. (4.61) 



zeCe 



fe=l 



The expansion (4.61) is uniform with respect to arg(^;) for \z\ — > oo in Ce and 

uniform in x as long as x varies in compact subsets of K intersected with the right 
Lebesgue set of B^^~^\ The expansion coefficients m+,fe(a;,ao) can be recursively 
computed from (4.5). 

Proof. To see that uniformity holds for this expansion, first recall the role of The- 
orem 3.11 in providing uniformity in the asymptotic expression (4.39) which then 
leads to (4.37) holding uniformly with respect to xq varying within compact subsets 
of M and with respect to arg (z) for |^!| — > oo in C^. This in turn leads to a similar 
uniformity holding for (4.46) which is the key to (4.57) holding with respect to xq 
varying within compact subsets of M and with respect to arg (z) for |z| — » oo in 

Ce. □ 

Remark 4.9. For simplicity, we focused thus far on the expansion of M+{z, xq, ao) 
as \z\ — > 00. Of course, Theorem 4.8 holds also for M-{z,xo,ao) replacing the 
hypotheses concerning right Lebesgue points by those of left Lebesgue points, etc. 
For convenience we just state the corresponding expansion and associated nonlinear 
recursion formula which covers both cases. 

N 

M±{z,x,ao) = ym±,k{x,ao)z-'' + o{\z\'^), N gN. (4.62) 



fe=o 



'm±fi{x,ao) = ±ilm, 

m±,i(x,ao) = -i(Bi,2(a;) +S2,i(a;)) ± ^(Bi^i^x) - B2,2{x)), 

k 



(K 
m'_j_ fc(a;, ao) + ^ m±,e{x, ao)m±,k+i-e{x, ao) 
e=i 

k 

+ X] '^±,^(2;, (XQ)B2,2{x)m±^k-i{x, ao) (4.63) 
1=0 

+ Si,2(a;)m±,fe(a;, ao) + m±,fe(a;, ao)B2,i(a;) ) , 
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l<k< N -1. 



Combining Theorem 4.8 and (2.77) then yields the analogous asymptotic expan- 
sion for M{z, X, ao). 

Theorem 4.10. Assume Hypothesis 2.1 with A = l2m, o^nd let ao = {Im 0) G 

(^mx2m p-^ e M and let a; e M. Suppose that for some N &n, B^^-^) e 
-^ioc(^)^'"^^™' X is a right and a left Lebesgue point of B^^~^\ and that for 
all R>0, 

rXo+R 



ess sup 

velxo,xo+R] 

+ ess sup 

ve[xo-R-xo] 

,= o{\zr 

\z\—>oo 



rXo+K -I 

/ dx' (a;') exp(2z^(:c' - y)) + —B^^-^^ (y) 

Jy 

r dx' S(^-i) (a;') exp(2i0(a;' - y)) - ^^B^^-^^ (y) 

Jx„-R 



C2 



(4.64) 



If N = 1, assume m addition Bk^k'Bgj, G LJ^JR)"^" for all k,k',l,l' G {1,2}. 
Let M{z,x,ao) be defined as in (2.76) {see also (2.77)). Then, as \z\ oo in C^, 
M(z..r, ao) has an asymptotic expansion of the form 

N 

M{z,x,ao) = {i/2)l2m+y2Mkix,ao)z-'' + o{\z\-''), N &N, 



(4.65) 



zeCe 



fe=i 



where 
Mi{x,ao) 



Bi,i{x + Q) - B2.2{x + 0) 

Sl,2(.T + 0) +B2,l(.T + 0) 
Bi,i(x-0)-B2,2(x-0) 

Si,2(a;-0) + B2,i(x-0) 



Bi,2(a: + 0)- 
S2,2(a; + 0) - 

Bi,2{x~0) 



•B2,i(.x + 0) 

Sl,l(.T + 0) 
-B2,l(x-0) 



52,2(a;-0)-Si,i(a;-0) 



etc. 
(4.66) 

The expansion (4.65) is uniform with respect to arg(z) for \z\ — > oo m and 
uniform, in x as long as x varies in compact subsets of R intersected with the right 
and left Lebesgue s et ofB^'^-'l 

If one merely assumes Hypothesis 2.1 with A = l2m, chq = {Im 0), and B G 

Lj„^(R)2'"x2'", then 



M{z,x,aQ) = {i/2)l2m + o{l). 

|z|— >oo 

zeCe 



(4.67) 



Again the asymptotic expansion (4.67) is uniform with respect to arg {z) for \z\ 
in Ce and uniform in x as long as x varies in compact intervals. 



oo 



The higher-order coefficients in (4.65) can be derived upon inserting (4.62) into 
(3.16a), taking into account (2.77). 

Theorems 4.7 and 4.8 (with AT e N) are new even in the scalar case m = 1 
with respect to the regularity assumptions on B. For previous results in the case 
m = 1 under stronger hypotheses on B we refer to [32], [56], [60], [61], [97]. In 
particular, [56], [60], and [61] derived alternative high-energy expansions for the 
Weyl-Titchmarsh m-fimction in the case m = 1. 

Throughout this section we fixed a to be ao = {Im 0). The case of general 
a e C^^x™ satisfying (2.9) then follows from (2.60). 
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B{x) 



(5.1) 



5. A Local Uniqueness Result 
In this section we assume that B is in the normal form given in Theorem 1.1, 

'Bi.i{x) Bi.2{x) 
,Si,2(a;) -Bi,i(a;), 

with Bi^\ and B\^2 self-adjoint a.e. We prove fundamental new local uniqueness re- 
sults for B in terms of exponentially small differences of Weyl-Titchmarsh matrices 
M+{z, X, a) and M(z, x, a). These results, in turn, yield new global ramifications. 
We start with an auxiliary result concerning asymptotic expansions. 

Lemma 5.1. Suppose a = {ax a^) G C™^^™ satisfies (2.9), fix Xo,yo € 11^ with 
yo > xq, and let x > xq. Assume A = l2m, B e L^([a;o, oo))^™^^™, supp(i?) C 
[^Ojyo]; with B in the normal form given in (5.1) a.e. on (xo^yo). Then, the follow- 
ing asymptotic expansions hold for Q{z,x,xo,a), ^{z,x,xo,a), and U+{z,x,xo,a) 
associated with (2.2a), 



9(^;,a;,a;o,a) 



^{z,x,xo,a) 



-i{al + ia2) 



— Q2 + ict 

\z\^oc 2 \-i{-a2 + ia 
zec+ 



) exp(-i^;(x - xq)) (l + o(l)) , x > xq, 

(5.2) 

^ *\ ] exp(-«z(x - xo))(l + o(l)), X > xo, 

(5.3) 



U+{z,x,xo,a) = 

\z\ — *oo 



i{al -ia*2) ) ®^P(^^(^ ~ ^o))(l + ^ ^ (5-^) 



Next, we introduce the abbreviation 

C = —Bi^2 ~ iBi^i, C* 
and suppose in addition that 



-B 



1,2 



ess sup 



J dx' B{x') exp{2iz{x' - y)) + ^B{y) 



zep+ 



along a ray p+ C C+, and that 

-^i,i> -Si,2> Bi^iBi^2, Bi^2Bi^i & L^{[xo,yo]y 

Then, 



(5.5) 



o{\z\-'), (5.6) 



(5.7) 



Q{z,x,xo,a) 



z6p+ 



al + 10.2 
-i{a\ + 10.2) 



i_ f C{x){al+ia*2) 
Az \iC{x){al + iaj) 

C{x'YC{x'){al+ia^) 
-iC{x')*C{x'){al+ial) 



;f.f \ ( ^ ( -(A. + ^"1 

^\Z,X,Xn,a) = - •/ * , • *x 

zep+ 



_ C{xonat-ia*2) 
4:Z \-iC{xo)*{al - ia2, 



er"^'=-''°\l + o{\z\-'^)), 

X > Xq, 

l_ f C{xoy{~a*2-ial) ^ 
Az \-iC{xQ)*{-al- ial) 



(5.8) 
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4z \iC{x){-a*+ial) 

^zJxo \-tCix'yC{x')i-a*2+tal)J J ^ ' 

X > xo, (5.9) 

whenever xq is a right Lebesgue point of B and x is a left Lebesgue point of B, and 

{C{xr -C{xor){al-ta*^) ' 



k ly (fJ;w)-tr-»l)) +»(N-')). ^ > 

(5.10) 

whenever x is a right Lebesgue point of B. 

Proof. Since xq and a are fixed throughout this proof, we will temporarily suppress 
these variables whenever possible to simplify notations. Introducing 

6(2:, x) = 2Q{z, x) exp(iz(a; — Xq)), (5-11) 

the Volterra integral equation for 9 (cf. (4.8)), 

Qt^ x)=(^*^ cos(2;(a; - a;o)) + sm{z{x - xq))\ 
^ ' ' \q.2 cos{z{x — Xq)) — Q.\ sin(2;(a; — a;o)) ) 

- [ dx' K{z,x,x')JB{x')Q{z,x'), (5.12) 

Jxo 

can be rewritten in terms of that of Q in the form 
al + 10.2 \ , I c^i ~ *ci;2 



[ dx'{R{x')exp{2iz{x - x')) + S{x'))e{z,x'), (5.13) 



I Xo 

where we abbreviated 



Using the elementary algebraic facts 

(5.15) 

for any a,b G C"^™, iterating (5.13) yields 

al+ic4 \ , [ "i - ^"2 \ g2«(x-xo) 
^-i(at + ia^)y \i{al +ia^)j 

+ ^(-2)-™ / da (J?(6)e''^^"-^^^ + 5(a)) X 

m=l -^^0 

X / ' rf6 (i?(6)e''^^^^-«^^ + 5(6)) • • • (5.16) 
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J Xo 

in 



J Xo \ 



i{al — 10.2) . 
, c/t ^( a*i+ia*2 \ 2iz(u-xo)\ 



Applying the Ricmann-Lebesgue lemma to (5.16) then proves (5.2) assuming B G 

L"'^([a;o, 00))^™^^™, only. Assuming also (5.6) and (5.7) one can compute the next 
term in the asymptotic expansion (5.2) and then obtains (5.8) using (5.16) and the 
finite-interval variant of (4.1), whenever xo is a right Lebesgue point of B and x is 
a left Lebesgue point of B. 

Exactly the same arguments apply to <&. Introducing 

$(z, x) = 2^{z, x) exp{iz{x - Xq)), (5-17) 
the Volterra integral equation for <[>, 



^, s _ /— a2 cos(2;(a; — a;o)) + al sin(i2;(x — Xq))^ 
^ ^ ' \ a* cos(i2:(a; — a;o)) + a2 sin(2;(a; — a;o)) ^ 

- f dx' K{z,x,x')JB{x')^j(z,x'), (5.18) 

•/ Xo 

can be rewritten in terms of that of $ in the form 

-if iR{x')exp{2iz{x-x'))+S{x'))^{z,x'). (5.19) 

^ Jxo 

Iterating (5.19), taking into account (5.15), yields 



-i{-a2+ial)j \i{-a2 - iQ.\) 
+ ^(-2)-™ / da (ii(a)e''"("-«^) + 5(a)) X 

m=l "'^o 

X / d(,2{R{^2)e^"^^'-^'^ +S{^2))... (5.20) 

Jxo 

■ ■ ■ f^""' dU-i (i?(^„x-i)e2-(«'"--«--) + S{U-,)) X 

J Xo 

X (_-i-t.,) 

Applying the Ricmann-Lebesgue lemma to (5.20) the proves (5.3) assuming B G 

L^{[xo,oo))'^"^^'^"\ only. Assuming also (5.6) and (5.7) one can compute the next 
term in the asymptotic expansion (5.3) and then obtains (5.9) using (5.20) and the 
finite-interval variant of (4.1), whenever xo is a right Lebesgue point of B and x is 
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a left Lebesgue point of B. 

Finally, we turn to U+{z,x). Introducing 

V+{z, x) = U+{z, x) exp{—iz{x — xq)), (5-21) 

the Volterra integral equation for U+, 

i(ai -zai)) / dx' K{z,x,x')JB{x')U+{z,x'), 

(5.22) 

can be rewritten in terms of that of V+ in the form 

■ ^ \ 1 pyo 

i{a*~-7i)) + 2 X dx'{R{x')+Six')eM'2iz{x'-x)))V+{z,x'). 

(5.23) 

Iterating (5.23), taking into account (5.15), yields 



2«(«2-ii) X 



~ / n,* — in,* \ ^ ry° ry° 

rVo ryo 

''I. -"4)) 

ryo fVo 
X d^sS{^s)e''"^^'-^'^--- d^2eR{^2e)x 

X y^^^ 5(6..0 (,(:i:r4)) e««-^-«^^). (5.24) 

Next, we take into account the different normalizations of U+ and U+. Using 
U+{z,xq) = [I,n M+{z,xoYY (cf-, (2.69) and ^{z,xo,XQ,ao) = hm), one readily 
verifies the relationship 

U+^l{z,x) =U+^i{z,x)u+^i{z,Xo)~^ , U+^2{Z,X) =Uj^^2{z,x)u+^i{z,Xo)~^. 

(5.25) 

Thus, applying the Riemann-Lebesgue lemma to (5.24) then proves (5.4) (in agree- 
ment with (4.26)), assuming B G L^{[xo, 00))^™^^™, only. Assuming also (5.6) and 
(5.7) one can compute the next term in the asymptotic expansion (5.4) and then 
obtains (5.10) using (5.24) and (4.1), whenever x is a right Lebesgue point of B. □ 

In the special case m = 1 (and for a = (1 0)), the expansion (5.10) was stated 
in [54]. 

Next, we note an elementary result concerning the boundary data independence 
of exponentially close Weyl-Titchmarsh matrices. 

Lemma 5.2. Fixx^ G K and suppose Aj = hm, Bj = B* e Li([a;o, Xo+i?])^'"''^™ 
for all R > 0. Denote by M+j{z,x,a), x > xo, the unique Weyl-Titchmarsh 
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matrices corresponding to the half-line Dirac-type operators D^j{a), j = 1,2, in 
(2.88). Fix an a G C'"^^'" satisfying (2.9) and assume that for all s > 0, 

\\M+,i{z,xo,a)-M+,2{z,xo,a)\\cm.m = o(e-2M-)(«--)) (5.26) 

along some ray p+ C C+. Then, for all a G c»"x2m gQUgjyi^ig (2.9) and for all 
£ > 0, 

||M+,i(^,aro,a)-M+,2(^,a;o,a)||c-x- = ^(e-^i-W^"-^)) (5.27) 
along the ray p+. 

Proof. Using (2.57) and (2.60) one estimates 

\\M+^i{z,xo,a) - M+_2{z,xo,a)\\c^^x^ 

= \\M+^i{z,xo,Q.)* - M+^2(2:,a;o,a)||c'">^'" 

< \\[aa* — M+^i{z,xo,ayaJa*]~^\\c^x^ x 
X \\M+^i{z,xo,a) - M+,2(z,a;o,Q!)||c".x™ x 
X \\[aa* + aJa* M+^2iz,xo,a)]-^\\c^x^, (5.28) 

since by (2.13) 

aJa*aa* + aa*aJa* = 0, aa*aa* — aJa*aJa* = Im- (5.29) 
Moreover, since 

[aa* — iaJa*][aa* — iaJa*]* = Im, (5.30) 

by (5.29), one infers (5.27) from (5.28) and M+j{z,Xo,Oi) = iIm + o{l) as \z\ — > oo, 
z e C+, j = 1,2 (cf. (3.1)). ' □ 

Our principal new local uniqueness result for Dirac-type operators in terms of 
Weyl-Titchmarsli matrices then reads as follows. 

Theorem 5.3. Fix xq € R and suppose Aj = hm, Bj G L'^{[xo,xo + i?])2'"x2m 
for all R > 0. Suppose also that Bj is in the normal form given in (5.1) a.e. on 
{xq,oo), j = 1,2. Denote by Mj^^(z,x,a), x > Xq, the unique Weyl-Titchmarsh 
matrices corresponding to the half-line Dirac-type operators D^j{a), j = 1,2, in 
(2.88). Then, 

if for some a > 0, Bi{x) = B2{x) for a.e. x G {xo,xo + a), (5.31) 
one obtains 

||Mi,+ (^,a;o,a)-M2,+ (^,a;o,a)||cmx„. = 0(e-2M^)») (5.32) 

\z\^oo 

along any ray p+ C C+ with < arg(z) < tt and for all a € £mx2m gafigfyifig 
(2.9). Conversely, fix an a G c.rn.x2m gatisfying (2.9) and if m > 1, assume in 
addition that Bj G L°°{[xo, xq + o])^™^^™, j = 1, 2. Moreover, suppose that for all 

£ > 0, 

||Mi,+(^,xo,ai)-M2,+(^,xo,di)||c™x™ = o(e-2imW(--)), £ = i,2, 

\z\-yoo 

(5.33) 
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along a ray p+,i C C+ with < a.Tg{z) < n/2 and along a ray C C+ with 
7r/2 < arg(^) < tt. Then 

Bi{x) = B2{x) for a.e. x G [xq, xo + a]. (5.34) 

Proof. Since (5.32) follows from Theorem 4.6 and Lemma 5.2, it suffices to focus 
on the proof of (5.34). Moreover, applying Theorem 4.6, we may without loss of 
generality assume for the rest of the proof that 

supp(Bj) C [a;o, xo + a], j = 1, 2. (5.35) 

In the following, we will adapt the principal ingredients of a recent proof of the local 
Borg-Marchenko uniqueness theorem for scalar Schrodinger operators (i.e., for m = 
1) by Bennewitz [13], to the current Dirac-type situation. First we recall that by 
Lemma 5.2, (5.33) holds along the rays p+j, j = 1, 2 for all a = {ai 02) € C^^^m 
satisfying (2.9). To simplify notations in the following we will again suppress xq 
and a whenever possible and hence abbreviate, Q{z,x,Xo,a), ^{z,x,xo,a), and 
t/j,+ (z, .X, .xo, a) by Q{z,x), ^{z,x), and Uj^^{z,x), respectively. Next, denoting in 
obvious notation by 

*^^(-) = ■ *^(-) = feS:::!) • ''^.^(-) - te::!:;:!) ■ 

j = l,2,x>xo, (5.36) 

the solutions associated with Bj, j = 1,2, which are defined in (2.14b) and (2.17), 
we introduce 

gj^k{z,x) = (l)j^k{z,x)uj^+^k{z,x)*, j,k e {1,2}, X > xq. (5.37) 

Using the asymptotic expansions (5.2)-(5.4) for Qj{z,x), ^j{z,x), and Uj^+{z,x), 
and the analogous ones for Qj{z,x)*, ^j{z,x)*, and Uj^+{z,x)* , one verifies for 
each fixed x > xq, 

9j,k{z,x) = {i/2)Im + o{l), j,fce{l,2}, (5.38) 

|z|— >oo 

zec+ 

assuming Bj e L^([a;o, + i?])2"'x2m for all i? > 0, j = 1, 2, only. Next, using the 
fact that for each fixed x> xq, 

<i)i,k{z,x)-^(j)2,k{z,x) = Im + o{l), k = l,2, (5.39) 

|z|— »oo 

zec+ 

{ui,+,k{z,xy)-^U2,+,k{z,x)* = Im + o{l), k = 1,2, (5.40) 

|2|— >00 

zec+ 

by (5.3), (5.4), one concludes 

^l,k{z, x)U2,+,j{z, X)* - Ui,+,fc(-2, X)4>2,k{z, x)* 

= ^i,k{z,x)62,k{z,x)* - ei^k{z,x)(l)2,k{z,x)* 
+ <i)i,k{z,x){M2,+ {z)-M^^+{z))(j)2,k{z,x)* = o(l), (5.41) 

\z\ — ^oo 

using (5.38), (5.40), and M2,+ {z)* = M2,+ {z). Combining hypothesis (5.33) and 
(5.3), one infers 

\Wk{z,x){M2,+{z) - Mi^+{z))(l>2,k{z,x)*\\ = o(l), x&{xo,xo + a) (5.42) 

\z\^OQ 
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along the rays ^ = 1,2. Thus, (5.41) impHes 
\\(l)-i^k{z,x)62,k{z,x)* - 6i^k{z,x)^2,k{z,x)*\\ = o(l), xG{xo,xo + a) (5.43) 

\z\—^oo 

along the rays p+,e, i = 1,2. The analogous estimate (5.43) holds along the complex 
conjugate rays p+/, £ = 1,2, in the lower complex half-plane C_. To simplify 
notations we denote the open sector generated by p+^i and its complex conjugate 
p+.i by Si, the open sector generated by the p+_2 and its complex conjugate p_|_.2 
by 52, the remaining sector in C+ is denoted by S3, and its complex conjugate 
sector in C_ is denoted by S4. Thus, one obtains a partition of C into 

4 

C=\JS'e, (5.44) 
t=i 

where each sector 5^, 1 < £ < 4, has opening angle strictly less than tt. Since 
(each matrix element of) the expression under the norm in (5.43) is entire and of 
order less or equal to one, one can apply the Phragmen-Lindelof principle (cf., e.g., 
[100, No. 322, p. 166 167, 379]) to each sector 5^, 1 < £ < 4, and obtains that 
each matrix element under the norm in (5.43) is uniformly bounded in each sector 
and hence on all of C. By Liouville's theorem, these matrix elements are all equal 
to certain constants. By the right-hand side of (5.43), these constants all vanish. 
Thus, we proved 

4>i,kiz,x)02,k{z,x)* = 9i^k{z,x)4>2,k{z,x)* for aU x G {xo,xo + a) (5.45) 
and hence 

(pi,k{z,xy^0i^k(z,x) ^ 92.k(z,x)*{(j)2,k{z,x)*)~'^ for all X e {xo,Xo+a). (5.46) 
Differentiating (j)j^k{z,x)~^0j^k{z,x), j,k = 1,2, with respect to x yields 
{(l)j^i{z,x)-'^9j^i{z,x))' 

= 4)j^i{z,x)~'^{{Bj)i^i{x) - z){(f)j^2{z,x)(t)j,i{z,x)~'^9j^i{z,x) - 9j^2{z,x)), (5.47) 

{(l)j,2{z,x)-'^9j,2{z,x)y 

= (pj^2{z, x)~^ {{Bj)i^i{x) + z){(l)j,i{z,x)(t)j,2{z,x)~'^9j,2{z,x) - 9j,i{z,x)). (5.48) 

Multiplying (5.47) by (j).j^i{z,x)* (4)j,i(z,x)*y^ and using (2.93), (2.95), and simi- 
larly, multiplying (5.48) by (j)j^2{z,x)*{(j)j^2{z,x)*)-'^ and using (2.94), (2.96) then 
yields 

(</),- i(z,x)-i^?,-i(^,a;))' = 0,,i(z,.T)-i((B,)ia(x)-z)(0,-i(z,.T)*)-\ (5.49) 
{ct>j,2{z,x)-^9j,2{z,x))'^<P,,2{z,x)-\{B,)i,i{x)+z){(t>j,2(.z,xy)-\ (5.50) 

In exactly the same way one derives 

{9^,^{z,xr{<t>i,r{z,xr)-')' 

= {9j,i{z,xr{4.j,i{z,xr)-'cpj,2{z,xr - 92iz,xr)i{Bj)i,i{x) - z){ct>j,i{z,xr)-' 

= cj>j,,{z,x)-\{Bj),,,{x) - z){cl>j,,{z,xr)-\ (5.51) 

{9j,2{z,xnci>j,2{z,xr)-'y 

= {9j,2iz,xr{4.j,2{z,xr)-'cPjAz,xr - e,-,i(^,a;)*)((B,)i,i(a;) +^)(<^,-,2(^,x)*)-i 
= cj>j,2{z,x)-\{Bj),,,{x) + z){cl>j,2{z,xr)-\ (5.52) 
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using (2.89)-(2.92). Thus, (5.46) implies 








0i,i(^,a;)*((Bi)i,i(a;) - = 


02,1 




"V2,i(^,a;), 
(5.53) 


0i,2(^,a;)*((Si)i,i(x) + z)~^(j)iAz,x) = 


02,2 


,xn{B2)i,i{x) + zy 


"^2,2(2,2;), 
(5.54) 






a;)*((S2)i,i(x)-z)- 


(5.55) 




^2,2 


xr{{B2)i,i{x) + zy 


"'^2,2(^,0;) 

(5.56) 



for a.e. x e (a;o, xq + a). Thus far we only used Bj £ L^{[xo, Xo + ii])2™x2m fgj- a,ll 
R>0, j = 1,2 and (5.35). 

In the special ease rn = 1. each of the equations (5.53) (5.56) allows for the 
completion of the proof of (5.34). Indeed, using the fact that 



(l>j^k{z,x) ^ (j)j^k{z,x), 9j.k{z,x) = 9j^k{z,x), j,fcG{l,2}. (5.57) 
and taking for instance (5.53), one infers for a.e. x e {xo,xo + a), that 



biAz,x)^ _ {Bi)iAx) - z 



Mz,xY {B2)i,i{x) - z' ^^-^^^ 

Since all zeros (and poles) of the left-hand side of (5.58) have even multiplicity, 
while aU zeros (and poles) of the right-hand side of (5.57) are simple, one concludes, 
assuming only that Bj e L^{[xq,xo + i?])^^^ for aU i? > 0, j = 1, 2, that 

(-81)1,1(2;) = (52)1,1 (a;) for a.e. .T e (a;o, a;o + a). (5.59) 

Thus for the case m = 1, we see by (5.53), and (5.54), (5.57), and (5.59), for 
a.e. X e {xo,xo + a), that 

<t>lk{z,x) = cl>lkiz,x), k = l,2. (5.60) 

Now, (2.92), (5.46), and (5.57) show, for a.e. x G {xo,xo + a), that 

0i,i(z, a;) 01, 2(2, a;) 



0i,i(2;,a;)0i,2(-2:,a;) 

= 02,i(-2,a;)02,2(2;,a;). (5.61) 



9i,i(z,a;) 6*1,2(2,0;) 
__ <i>2.\{z,x) ^ 02, 2(2, a;) 
6*2, 1(2, a;) 6'2,2(-z,a;) 

By (2.2a) we see that 

(02i(z,x))' = 2(z- (B,(x))i,i)0,,i(z,a;)0,-2(z,a;) + (S,(a;))i,20,',i(z, x), 3 = 1,2. 

(5.62) 

Thus, by (5.59), (5.60), and (5.61), 

(Bi(a;))i,2 = (B2(a;))i,2 for a.e. x e {xo,xo + a). (5.63) 

Together, (5.59) and (5.63) imply (5.34) in the special case m = 1. 

Unfortunately, the case m > 1 appears to be quite a bit more involved. To deal 
with this case we first note that taking determinants in (5.53) yields 

det(0i,i(^,a;,a;o,Q!)*)det(0i,i(2;,a;,a;o,Q!)) det((Bi)i,i(a;) - zlm)) 



det{(t)2,i{z-:X,XQ,a)*)det{(j)2,i{z,x,XQ,a) det((B2)i,i(a;) - zlm)) 



(5.64) 



44 



CLARK AND GESZTESY 



for a.e. x € {xq, xq +a). Next, we intend to prove that 

det((Bi)i,i(a;) - zim)) = det((B2)i,i(a;) - zim)) for a.e. x e {xo,xo + a). (5.65) 

Given the fact that {Bj)i^i{x), j = 1, 2, is self-adjoint, showing (5.65) is equivalent 
to showing that Bi{x) and B2{x) are unitarily equivalent for a.e. x € {x(j,XQ + a). 
Arguing by contradiction, we assume that at least one pair of eigenvalues of Bi {x) 
and B2{x) differs. Thus, fixing xi € {xo,xo + a), let A(.ti) be an eigenvalue of 
Bi{xi) but not ofB2{xi). Then (5.64) implies, for aU a € C^^^m satisfying (2.9), 
that 

det(0i,i(A(a;i),a;i,a;o,Q;)) =0. (5.66) 
Next, for A e M and x > xq define 

N{X, X, a) =^i,2(A, X, xo,a)di^2{\ x, xq, a)* 

+ 0i,2(A,a;,a;o,Q!)0i,2(A,a;,a;o,a)*. (5.67) 
Then, A''(A, x, a) is strictly positive definite, 

A^(A,x,a)>0. (5.68) 
Indeed, suppose Nf — for some / G C™, then 

^i,2(A)ei,2(A)*/ + 0i,2(A)0i,2(A)*/ = (5.69) 



implies 
and hence 
Thus, 



9i,2(A)0i,2(A)*/ = 0, (/.i,2(A)(/.i,2(A)*/ = (5.70) 
^i,2(A)*/ = 0, <^i,2(A)*/ = 0. (5.71) 



/ = (^i,i(A)<^i,2(A)* - MX)6i,2{Xy)f = (5.72) 

by (2.96), and hence / = proves (5.68). Introducing ao = (/„ 0) e C^^^™ and 

7 = (71 72) G C"^2m defined by 

71 = [Oi.2{H^i),xi,XQ,aQ)9i,2{X{xi),xi,XQ,aQ)* (5.73) 

+ (^i,2(A(a;i),a;i,a;o,ao)</>i,2(A(a;i),a;i,a;o,Q;o)*]"^/^6'i,2(A(a:i),a;i,a;o,ao), 

72 = [6'i,2(A(a:i),a;i,a;o,ao)6'i,2(A(a;i),a;i,a;o,ao)* (5-74) 

+ 4>i^2{X{xi), xi,xo, ao)4>i,2iX{xi),xi,xo, aoy]~^^'^4>i,2{X{xi), xi,xo, ao), 

one verifies 77* = 7„ (by (5.73) and (5.74)) and 7J7* = (by (2.94)). Thus, 7 
satisfies (2.9). Next, since 

<?^i,i(A(a::i),a;i, 2:0,7) = (/)i,i(A(a:i),Xi,Xo,Q!o)7i - ^i,i(A(a:i), Xi, xq, 00)73 (5.75) 
as a special case of (2.97), one derives 

(l)i,i{X{xi),xi,xo, 7) = [(/)i,i(A(xi), xi,xo, ao)6'i,2(A(a;i), xi,xo, ao)* 

- 6'i,i(A(a:;i),a;i,a;o,ao)0i,2(A(a;i),a:i,a;o,Q;o)*] x 
X [0i^2{X{xi),xi,xo,ao)9i,2{X{xi),xi,xo, ao)* 

+ 4)i^2{X{xi), xi, xo, Q;o)<?i'i,2(A(a;i), xi,xo, ao)*]~^^^ 
= -[^i,2(A(a;i), xi,xo, ao)0i^2{X{xi), xi,xo, ao)* (5.76) 
+ (?!>i,2(A(a;i),a;i,a;o,ao)(Ai,2(A(a;i),a;i,a;o,Q!o)*]~^/^ < 0. 
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using (2.96). This contradiction to (5.66) proves (5.65). Hence for A e M and for 
a.e. X e {xo,xo + a) 

|det((^i,i(A,a;,a;o,a))| = \det{(l)2,i{X,x,xo,a))\, (5.77) 

by (5.64). Equation (5.77) implies that for a.c. xi G (a;o,a;o + a), the family of Dirac 
operators £'+(a, ao) in L^([a;o, Xi])^™, defined by 

D+{a,ao) = J^-B, (5.78) 

dom(£>+(a,ao)) = {<A e L\[xo,xi]f"' | <p e AC([a;o,a;i])2'"; 

a<f>{xo) = 0, ao<l>ixi) = 0; (J<^' - B(f>) e L''{[xo,xi]f"'}, 

with ao = {Irn 0), have identical spectra for all boundary data a E C™x2m g.^^jg_ 
tying (2.9). Hence, assuming Bj S L°°{[xo,xo + a])^™^^™, j = 1,2, one can apply 
Theorem 2.3 of Malamud [88] and obtains (5.34). □ 

We should note that Malamud's Theorem 2.3 in [88] only requires the equality 
of + 1 spectra (associated with linearly independent boundary data indexed by 
a G C™^2m) in to conclude (5.34). 

There is no particular significance of the rays pi, £ = 1,2, in Theorem 5.3. Any 
non-selfintersecting Jordan arc that tends to infinity in the sectors e < aTg{z) < 
(7r/2) — £ and (7r/2) + e < arg(^) < tt — e for some < £ < 7r/4 will do. 

Remark 5.4. We were not able to prove (5.34) directly from (5.53)-(5.56), with- 
out resorting to the arguments involving (5.77) and (5.78). To conclude the proof 
according to the Borg-type Theorem 2.3 of Malamud [88] (cf. also Theorem 4 in 
[89]), requires the introduction of the extra hypothesis Bj G L°°{[xo, .To + a])^™^^™, 
j = 1, 2 in the matrix context m > 1, since the construction of transformation op- 
erators for Dirac-type systems, to date, uses such an additional hypothesis on B. 
This extra hypothesis is clearly superfluous in the case m = 1 . Obviously, one con- 
jectures that this extra hypothesis on Bj should also be redundant in Theorem 5.3, 
but this appears to require nontrivial future efforts. In this context it might be 
interesting to note that the higher-order expansions (5.8)-(5.10) do not determine 
B imiqucly. An explicit analysis shows that while they do determine -Bi,2) they 
only determine B^ i, not itself. So that approach does not aide in proving 
(5.34) (besides, it would require the additional hypotheses (5.6) on B). 

The corresponding local uniqueness result in terms of M{z, xq, a) then reads as 

follows. 

Theorem 5.5. Fix xq € R and suppose Aj = hm, Bj € Ll^^{U.)'^'^^'^"' , and 
Bj = Bj a.e. on R, j = 1,2. Suppose also that Bj is in the normal form given in 
(5.1) a.e. on {xq,oo), j = 1,2. Denote by Mj{z,xa,a), the unique Weyl-Titchmarsh 
matrices (2.76) corresponding to the Dirac-type operators Dj, j = 1,2, in (2.75). 
Then, 

if for some a > 0, Bi{x) = B2{x) for a.e. x G {xq — a,xo + a), (5.79) 
one obtains 

\\Mi{z,xo,a)-M2{z,xo,a)\yr..2^ = ©(e'^M^)") (5.80) 
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along any ray p+ C C+ with < a.vg{z) < n and for all a G C"*^^"* satisfying (2.9). 

Conversely, fix a a G C™^^™ satisfying (2.9) and if m > 1, assum,e in addition that 
Bj e L°°{[xo - a,xo + a])2"^2", j = 1, 2. Moreover, suppose that for all e > 0, 

||Mi(^,a;o,ai)-M2(z,xo,di)||c-x2„ = o(e-2M.)(a-e)^^ £ = 1,2, (5.81) 

a/ong' a ray p+^i C C+ with < arg(2;) < n/2 and along a ray C C+ with 
7r/2 < arg(2;) < tt. T/ien 

Bi{x) = B2{x) for a.e. x G [xq — a,Xo + a]. (5.82) 

Proof (5.80) is proved by combining (2.77), and (5.31), (5.32), and (5.82) then 
follows by combining (2.77), and (5.33), (5.34), taking into account the asymptotic 
expansions 

M±{z,xq) = ±ilm + 0{l) (5.83) 

|z|— >oo 

along any ray with £ < arg(^;) < tt — e in the case of Dirac-type operators (cf. (3.1)). 

□ 

Remark 5.6. Theorem 5.3 and Theorem 5.5 yield new global imiqucness theo- 
rems for half-line and full-line Dirac-type operators, extending the classical Borg- 
Marchenko-type results. Indeed, if (5.33) (resp., (5.81)) holds for all a > 0, then 
(5.34) (resp. (5.82)) holds for a.e. x G [a;o,oo) (resp., for a.e. x e M). 

In the case of scalar Schrodinger operators, the analog of Theorem 5.3 is due to 
Simon [114]. An alternative proof, applicable to matrix- valued Schrodinger opera- 
tors was presented in [50] (cf. also [41]). More recently, yet another proof was found 
by Bennewitz [13] (following some ideas in [16]). These results extend the classi- 
cal (global) uniqueness results due to Borg [16] and Marchenko [91], [92] (cf. also 
[14]), which state that half- line m- functions uniquely determine the corresponding 
potential coefficient. The Dirac-type results presented in this section (especially, 
all local considerations) appear to be new, even in the special case m = 1. Pre- 
vious results in the Dirac case focused on global uniqueness questions only. We 
refer to Gasymov and Levitan [34] in the case m = 1 and to Lesch and Malamud 
[81] in the matrix case m G N. Most recently, Alexander Sakhnovich kindly in- 
formed us that his integral representation of the Weyl-Titchmarsh matrix in [103] 
can be used to derive asymptotic expansions for the Weyl-Titchmarsh matrix and 
its associated matrix- valued spectral function, and also yields a result analogous to 
Theorem 5.3 (i) for a certain class of canonical systems. Moreover, in the case of 
skew-adjoint Dirac-type systems, similar results are discussed in [104] and applied 
to the nonlinear Schrodinger equation on a half-axis. 

Although not directly used in this paper, it should be pointed out that inverse 
monodromy problems for canonical systems received a lot of attention (some of it 
very recently). The interested reader is referred to [4], [5], [6], [87], [88], [89], [109], 
[112] and the extensive literature cited therein. Moreover, inverse spectral theory 
associated with canonical systems is discussed in [96], [104], [106], [107], [109], [110], 
[111], [112] (see also the extensive literature cited in [41]). 

6. Trace Formulas and Borg-Type Theorems 

In our final section we derive a trace formula for B and then discuss its applica- 
tion to Borg-type uniqueness theorems for Dirac-type operators. 
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Theorem 6.1. Assume Hypothesis 2.1 with A = l2m, and let ao = {Im 0) e 
Cmx2m_ p^j. xq^^ and suppose that for all R>0, 



ess sup 

ye[xo,xo+-R] 



CSS sup 
ye[xo-R,xo] 



dx' B{x') exp(2i^(a;' - y)) + ^^B{y) 
r dx' B{x') exp(2i0(a;' - y)) - ^B{y) 



C2mx2n 



C2mX2n 



|z|— >oo 



oi\zn 



along a ray p+ C C+. In addition, assume Bk.k'Bt.i' 
k,k',£,£' S {1,2}. Then, with T(A,x, ao) defined in (2.83), 

'Si,i(x) - B2,2{X) Bi,2{x) + B2,l{x) 

^B^,2{x) + B2,i{x) B2,2{x) - Bi,i(a;) 



(6.1) 
for all 



lim 2 / dXz^iX- z)-^T{\,x,ao) f 



1^ 



or a.e. x e 



(6.2) 



Proof By (2.78), 



^ln{M{z,x,ao)) = [ d\{\ - z)-^r{\,x,ao). (6.3) 

dz 

Next, suppose that a; e M is a left and right Lebesgue point of B. By (4.65), (4.66) 
one obtains 
d 



dz 



]ii{M{z,x,ao)) 



1 /Si,i(a; + 0)-B2,2(a; + 0) Bi,2(a; + 0) + B2,i(a; + 0) 
I^Hoo 4 \Bi^2{x + 0) + B2,i(a; + 0) ^2,2(0; + 0) - Bi,i(a; + 0) 



2Gp+ 



1 - 0) - B2,2{x - 0) Bi,2(a; - 0) + ^2,1(2; - 0) 



4 \Bi,2{x - 0) + S2,i(a; - 0) B2,2(a; - 0) - - 0) 



(6.4) 

+ o(z-) 



and hence 



1 + 0) - ^2,2(2; + 0) Si,2(x + 0) + 52,i(a: + 0) 

2 V^i,2(a; + 0) + S2,i(a; + 0) 52,2(2; + 0) - Bi,i(a; + 0) 

1 /Si,i(,x-0) -53.2(3; -0) Bi,2(x-0) + B2.i(a;-0) 

2 \Bi,2{x - 0) + £2,1(2; - 0) £2,2(2; - 0) - Bi,i(a; - 0) 



= Hm 2 / dA^2(;^_^-)-2y(;^^^^^^)_ 



Sincc a.e. .t G M is a Lobcsguc point of B, one concludes (6.2) 



(6.5) 



□ 



In the case to = 1, a trace formula for Dirac-type operators, using Krein spectral 
shift functions and exponential representations of Herglotz functions, was discussed 
in [116]. This circle of ideas was first introduced in connection with trace formulas 
of Schrodinger operators in [48] (see also [38], [39], [101], [102] in the scalar case 
TO = 1. The corresponding case of trace formulas for matrix- valued Schrodinger 
operators was introduced in [37] (sec also [21]). 

Analogous trace formulas can be drived for all higher-order coefficients Mk{x, ao) 
in (4.65) (see, e.g., [39] in connection with scalar Schrodinger operators). 
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A comparison of the trace formula (3.20) in [21] for Schrodinger operators with 

its Dirac-type counterpart (6.2) reveals characteristic differences. While in the 
Schrodinger case the trace formula directly involves the potential coefficient Q{x), 
Mi{x,ao) differs markedly from a constant multiple of B{x), and consequently, 
the Dirac-type trace formula (6.2) does not directly involve B{x) but certain linear 
combinations of Bj^k{x). This is related to the fact that M[z,XQ,ao) (or equiva- 
lently, T(A, ccq, ao))) in general, does not uniquely determine B a.e. In fact, there 
exists a typical ambiguity concerning the coefficients of D related to unitary gauge- 
transformations oiD. In the case m = 1 this ambiguity is well-known and discussed, 
e.g., in [34], [84, Sect. 1. 10], [85, Ch. 7]. These gauge transformations leave the spec- 
trum of D invariant and suggest that we focus our attention on certain normal forms 
of D in connection with inverse spectral problems for Dirac-type operators. 

Lemma 6.2. Assume Hypothesis 2.17. Then D = J-^ — B is unitarily equivalent 
to D, where D in L^(R)^'" is of the normal form 

D = J±-B={ /^'l -I^±-B.,.\ 

Here B = B* a.e. and 

= -(l/2)Im(;7i;^[(Si.2 + B2,i) - i(Bi,i - B2,2)]C/2,2) = i, (6.7) 

i?i,2 = il/2)Re{U{-j[{Bi,2 + B2.1) - x(Bi,i - B2,2)]C/2,2) = i?i*,2, (6.8) 

with Ujj € ([ynxm ^ j = 1,2, satisfying the first-order system of ordinary differential 
equations 

iU'^jix) = -{l/2){{-iy{B,,i{x)+B2,2{x))+i{BM - B2,,{x)))Ujj{x), 

for a.e. x€R, j = 1,2. (6.9) 

Proof. We start with the unitary transformation V in i^(]R)^'" defined by 



v = ^[ 7" .7 , v-' = ^i ^7 , (6.10) 



which maps D to Di, where 

d 



Di = V-'DV 



^dx " ^ 



_ 1 / Bi^i + i?2,2 — *(-Bl,2 — -62,1) i?l,2 + -62,1 — i{B\^i — i?2,2) 

2 \B^,2 + -82,1 + i(Bi,i - -62,2) -Bi,i + ^2,2 + i(-Bi,2 - -B2,i) 
Next, we introduce the unitary operator U in i^(M)^'" defined by 



(6.11) 



U=i^l;' ° J, (6.12) 



U2 



where the unitary m x m matrices Uj,j E c"ix™ ^re solutions of the first-order 
system (6.9). Since by hypothesis Bj^k e Lj„^(R)'"^™, l<j,k<2, the solutions of 
equation (6.9) are well-defined and Ujj G ACioc(R)'"^'", j = 1,2. One computes 

D = U-^D,U=(''-f -^''l), (6.13) 

\~^i,2 ""was/ 
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where Bi,2 e Lf„^(M)'"^™ and 

Bi,2(a;) = {l/2)U^l{x){Bi^2{x)+B2,i{x)-iiBi^i{x)-B2,2{x))U2,2ix)- (6-14) 

Finally, defining D = VDV~^, one arrives at (6.6)-(6.8). □ 

Thus, unitary invariants of D (such as the spectrum, spec(£>), of D and its 
multiplicity) cannot determine B in general but at best a potential matrix of the 

type (normal form) B in (6.6). A further restriction on the solvability of inverse 
spectral problems for Dirac-type operators is mentioned in the following result. 

Lemma 6.3. Assume Hypotheses 2.17 and let co = lu* £ C™^™ be a constant self- 
adjoint mxm matrix. Then D = J-^ — B is unitarily equivalent to inL^ (R) ^"^ , 
where 

- Ba,,l,2 J 

with 

^^,1,1 = -(l/2)Im(e'"C/i:iM(Si,2 + B2,i) - - S2,2)]!72,2e*'") = 

5^,1,2 = {l/2)Rc{e'^U^j[iBi.2 + B2,i) - i{Bx^i - 52,2)1(72,26*'") = S^^^^j, (6-16) 

and with Ujj, j = 1,2, satisfying the first-order system (6.9). 

Proof. Define 

Uo.= (^7 AV (6.17) 




^ 

Using the notation employed in the proof of Lemma 6.2 one verifies that 

= VU^{VU)-^DVU{VU^)-\ (6.18) 

□ 

In particular, choosing ui = {n/2)Im effects the sign change B —B, with B 
given by (6.7), (6.8). 

For detailed discussions of various normal forms for Dirac-type operators we refer 
to [34], [59], [84, Ch. 9], [85, Ch. 7] in the case m = 1 and to [33], [81], [88], [93, 
p. 193-195], [95] in the general matrix- valued case. Perhaps it should be noted that 
if D is in its normal form D as in (6.6), turns into a 2m x 2m matrix- valued 
Schrodinger operator under appropriate regularity assumptions on B. Details on 
this fact and the relation between the M-matrices of D and D"^ can be found in 
Section 3 of [41]. 

Next, we turn to Borg-type theorems, one of the principal topics of this paper. 
In 1946 Borg [15] proved, among a variety of other inverse spectral theorems, the 
following result. 

Theorem 6.4 ([15]). Assume q S Ll^^{R) to he real-valued and periodic and let 

h = -d'^/dx'^ + q (6.19) 

be the associated self-adjoint Schrodinger operator in L^(R). Moreover, suppose 
that spec(/i) = [eo, oo) for some bq € K. Then 

q[x) = 60 for a.e. x€ M. (6.20) 
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The analog of Theorem 6.4 for Dirac-type operators (in the case m = 1) was 

proven by Giacheti and Johnson [53] in 1984 (sec also [35], [36], [47] in the special 
case where p is constant and [55] in the case where p,q G L"^ (M) are real- valued and 
periodic). 

Theorem 6.5 ([53]). Assume p.q ^ L°°{R) to be real-valued and periodic and let 




(6.21) 



be the associated self-adjoint Dirac-type operator in L^(M)^. Moreover, suppose that 
spec(d) = R. Then 

p{x) = q{x) = for a.e. xe M. (6.22) 

Traditionally, uniqueness results such as Theorems 6.4 and 6.5 are called Borg- 
type theorems. (However, this terminology is not uniquely adopted and hence 
a bit unfortunate. Indeed, inverse spectral results on finite intervals recovering 

the potential coefficient (s) from several spectra, were also pioneered by Borg in 
his celebrated paper [15], and hence are also coined Borg- type theorems in the 
literature, see, e.g., [86], [88], [89].) 

A quick and natural proof of Theorem 6.4, based on a trace formula for q, was pre- 
sented in [21]. This strategy of proof was then applied to the case of matrix- valued 
Schrodingcr operators and the corresponding matrix-valued analog of Theorem 6.4 
was also proved in [21] along these lines. A closer examination of the proof of The- 
orem 6.4 shows that periodicity of q is not the crucial element in the proof of the 
uniqueness result (6.20). The key ingredient (besides spec(/i) — [eo,oo)) is clearly 
the fact that for all a; e K, 

^{X,x) = 1/2 for a.e. A e specess(/i) (6.23) 

(specess( • ) the essential spectrum), where ^{-jx) is defined by 

^(A, x) = lim 7r"^Im(ln(5f(A -|- ie, x))) for a.e. A e R, (6.24) 

£—►0 

and g{z,x) denotes Green's function (i.e., the integral kernel of the resolvent) of h 
on the diagonal, 

g{z,x) = {h- z)~'^{x,x). (6.25) 

Completely analogous considerations apply to the Dirac-type case. 

Real- valued periodic potentials are known to satisfy (6.23) but so arc certain 
classes of real- valued quasi-periodic and almost-periodic potentials q (see, e.g., [23], 
[24], [25], [70], [72], [75], [76], [77], [78], [115]). In particular, the class of real- valued 
algcbro-gcomctric finite-gap potentials q (a subclass of the set of real- valued quasi- 
periodic potentials) is a prime example satisfying (6.23) without necessarily being 
periodic. Traditionally, potentials q satisfying (6.23) are called reflectionless (see 
[24], [25], [77], [115]). Again the analogous notions apply to the Dirac-type case 
(cf., e.g., [23], [53], [71]). 

Taking this circle of ideas as the point of departure for our derivation of Borg- type 
results for Dirac-type operators, we now use the reflectionless situation described 
in (6.23), actually, its proper analog for Dirac-type systems, as the model for the 
subsequent definition. 
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Definition 6.6. Assume Hypothesis 2.1 with A = l2m, and let ao = {Im 0) G 
Cmx2m rpi^gjj 5 jg p^ugjj reflectionless if for all xeR, 

T{X,x,ao) = (l/2)72m for a.e. A e speCess(-D). (6.26) 

Since hardly any confusion can arise, we will also call the Dirac-type operator D 
reflectionless if (6.26) is satisfied. 

Given Definition 6.6, we turn to a Borg-type uniqueness theorem and formulate 
the analog of Theorem 6.4 for (reflectionless) Dirac-type operators. 

Theorem 6.7. Assume Hypothesis 2.1 with A = l2m, and let ao = (Im 0) G 
(^mx2m_ jj: X G R, T{X,x,ao) = C is a constant 2m x 2m matrix for a.e. 

X gR, especially, if B is reflectionless and spec(D) = R, then 

Bi,i(a;) = B2,2{x), Bi^2{x) = -B2,i{x) for a.e. x &R. (6.27) 

In particular, if D is assumed to be in its normal form (6.6), that is, of the type 
D = jA.-B, then 

B{x) = for a.e. xG R. (6.28) 

Proof The fact that J^dX (X - z)'^ = for all z e C\M, that a.e. a; e M is a 
Lebesgue point of B, and the trace formula (6.2), imply (6.27). Together with 
Lemma 6.2 this yields (6.28). □ 

The analog of Theorem 6.7 for matrix- valued Schrodinger operators was recently 
proved in [21]. 

In the remainder of the section we will show that the case of periodic B is covered 
by Theorem 6.7 under appropriate uniform multiplicity assumptions on spec(D). 
In order to handle Floquet theoretic aspects of periodic Dirac-type operators D, 
we adopt the following assumptions until the end of this section. 

Hypothesis 6.8. In addition to Hypothesis 2.1 assume A ~ l2m and suppose that 
B is periodic, that is, there is an w > such that B(x + uj) = B{x) for a.e. x & R. 

The following result has been proven in [21, Theorem 4.6]. 

Theorem 6.9 ([21], Theorem 4.6). Assume Hypothesis 6.8 and let = (/„ 0) G 
^2mxm_ Jj JJ j^^g uniform spectral multiplicity 2m, then for all x G R and all 
X e spec(£>)°, 

M+{X + iO,x,ao) = M_{X + iO,x,ao)* = M_{X - iO,x,ao). (6.29) 

In particular, M-{z,x,ao) is the analytic continuation of M+{z,x,ao) {and vice 
versa) through spec{D)° . 

Here A° denotes the open interior of a set C M. 

Strictly speaking. Theorem 4.6 in [21] was proved for matrix-valued Schrodinger 
operators. But the proof extends line by line to the corresponding Dirac-type 
situation and was predominantly formulated in terms of Hamiltonian systems no- 
tation (rather than Schrodinger operator specifics) in order to be applicable to the 
present context. In particular, the spectrum, spec{H), of the Schrodinger operator 
H should be replaced by that of D, the point spectrum, speCp(iJ^), of the Dirich- 
let Schrodinger operator H^^ with a Dirichlet boundary condition at the point xq 
should simply be replaced by the set {A e M | det((/ii(A, xq + oj, xo, ao)) = 0}, etc. 
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Theorem 6.10. Suppose Hypothesis 6.8 and let ao = {Im 0) e C^™^™. If D has 
uniform spectral multiplicity 2m, then D is reflectionless and for all x gM. and all 
A G spec(Li)°, 

T(A,x,ao) = (l/2)/2m. (6.30) 
Proof. This is clear from (2.77) and (6.29), which imply 

M(A + iO, X, ao) = -M(A + iO, x, ao)*. (6.31) 

□ 

Theorems 6.9 and 6.10 extend to more general situations (not necessarily periodic 
ones) as is clear from the corresponding results in [23], [53], [42], [75], [76], [77], 
[115] in the scalar case m = 1 (replacing the phrase "for all A G spec(£')°" by 
"for a.e. A e spec(£')°", etc.). For the corresponding matrix-valued Schrodinger 
operator case we refer to [78]. 

Corollary 6.11. Assume Hypothesis 6.8. If D has uniform spectral multiplicity 
2m and spec{D) — R, then 

Bi,i{x) = B2.2{x), B,Jx) = -B2,i(x) for a.e. x e R. (6.32) 

In particular, if D is assumed to be in its normal form D = J-^ — B, with B given 
by (6.6), then 

B(x) = for a.e. a; G M. (6.33) 

Remark 6.12. The assumption of uniform (maximal) spectral multiplicity 2m in 
Corollary 6.11 is an essential one. Otherwise, one can easily construct nonconstant 
potentials B such that the associated operator D has overlapping band spectra 
and hence spectrum the whole real line. Also self-adjointness of B is crucial for 
Corollary 6.11 to hold (cf. the corresponding discussion in Remark 4.2 of [21] in the 
context of Schrodinger operators) . 

The analog of Corollary 6.11 for periodic matrix- valued Schrodinger operators 
was first proved by Dcprcs [26] and recently rederived using such a trace formula 
approach in [21]. 

We note that all results presented in this paper also apply to matrix-valued 
finite-difference Hamiltonian systems. We refer the reader to [22] in this direction. 

Finally, Borg-type uniqueness theorems for Hamiltonian systems are just a be- 
ginning. There is a natural extension of Borg's Theorem 6.4 to self-adjoint periodic 
Schrodinger, respectively, Dirac-type operators with one gap, respectively, two gaps 
in their spectrum. In the case of (scalar) Schrodinger operators, such an extension 
is due to Hochstadt [69] and the resulting potential q becomes twice the elliptic 
Weierstrass function. In the case of Dirac-type operators (with m = 1 and vanish- 
ing diagonal coefficients in B) such an extension involving elliptic functions can be 
found in [35], [36], [47] (see also [52]). Extensions to matrix- valued versions (i.e., 
for m > 2) are currently under active investigations. 
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